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Zut Theorie der Mahlvorginge 


Von Lupwie Bass, Reutte, Tirol?) 


Einfihrung 


Die grundlegende Frage der vorliegenden Theorie der Mahlvorginge besteht 


_ darin, bei gegebenen, festgehaltenen Mahlbedingungen die Gewichtsverteilung 
_ des Mahlgutes als Funktion der Mahldauer auszudriicken. Es soll also bei be- 
_ kannter Anfangsverteilung M(x, 0) die Funktion M(x, ¢) bestimmt werden, 
_ die dadurch definiert ist, dass der Ausdruck 


M(x, t) dx 


| das Gewicht der Korner des Mahlgutes, deren Gréssen zwischen x und x + dx 
 liegen, zur Zeit ¢ angibt. 


Die vorliegende Theorie ist im wesentlichen eine Prat omen slaseiia: nur 
so konnten die. komplizierten und experimentell wenig erforschten Vorgiénge 
in Mihlen verhaltnismassig einfach mathematisch erfasst werden. Die Theorie 
beruht insofern auf statistischer Grundlage, als die dort eingefiithrten Gréssen 
erst bei einer grossen Anzahl der Korner aller im Mahlgut vorhandenen Grés- 


_ senbereiche reproduzierbar werden. 


Es sei e(€, x) die Masse der Kérner im Gréssenintervall (x, x + dx), die beim 


_ Mahlen pro Zeiteinheit aus einer Einheitsmasse der Kérner des Grdésseninter- 
_valls (¢, € + d&) entstehen. Dieser Definition treten drei Voraussetzungen bei. 


1. Bei dem Mahlen geht eine «reine Zerkleinerung» vor sich, das heisst, dass 


‘sich kleinere Kérner nicht zu grésseren vereinigen. Daraus folgt, dass e(€, x) 


nur im Bereich 0 <x S&<x,, von Null verschieden sein kann (mit x, 
wird die maximale im Mahlgut vorhandene Korngrosse bezeichnet). 
2. Die Funktion «(&, x) ist zeitunabhangig. Diese Voraussetzung ist nur 


dann berechtigt, wenn die durch die Korngréssenveranderungen verursachten 


Effekte, vor allem das Polstern der Miihle mit sehr kleinen K6rnern, vernach- 


-_ladssigt werden kénnen (das Volumen des Mahlgutes muss sehr klein sein gegen- 


iiber dem Volumen der Miihle). 


1) Metallwerk Plansee. 
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3. Die Funktion e(&, x) ist im Bereich 0 S % S — < x,, stetig. Die Berech- 
tigung dieser Annahme beruht im folgenden: Es kénnen zwar plétzliche starke 
Verdinderungen von e(€, x) in solchen Bereichen erscheinen, wo sich die Eigen- 
schaften des Materials und der Mahlvorrichtung, die das Zerspalten der K6rner 
beeinflussen, plétzlich andern. Ein solches Verhalten von e(&, x) in der Koor- 
dinate x ist zu erwarten, wenn die Kérner bis zur Grésse der sogenannten 
idealen Gitterbereiche!) zerkleinert werden, bei der ihre Festigkeit stark zu- 
nimmt, und ebenfalls in der Koordinate &, wenn die Kérner die Grésse der 
Poren der Mahlkérper erreichen. Nun besteht das Mahlgut aus sehr vielen 
idealen Gitterbereichen, deren Grdésse in einem gewissen Intervall schwankt; 
die Verteilungskurve dieser Gitterbereiche wird deshalb annahernd stetig, so 
dass, auch wenn die Idealkristalle in der Miihle nicht zerspaltet werden kénnen, 
e(é, x) stetig bleibt. Die zahlreichen Poren der Mahlkérper sind ebenfalls nicht 
alle gleich gross, so dass sie bloss um so starkere Veranderungen von é(, x) in 
€ bewirken kénnen, je kleiner ihr Gréssenbereich ist. 

Ahnliche Uberlegungen fiihren zur Annahme der Stetigkeit von M(x, #). 
Da M(x, t) im ganzen Korngréssenbereich definiert und endlich ist, sind nur 
Unstetigkeiten durch endliche Spriinge denkbar. Zu solchen ware aber eine 
Anzahl von Kérnern genau gleicher Grésse notwendig; dies trifft aber nicht zu, 
da sich praktisch jedes Korn im Gewicht vom andern etwas unterscheidet und 
da die Theorie eine grosse Anzahl von Koérnern aller vorhandenen Gréssenbe- 
reiche voraussetzt. Anstatt der Sprungstellen kénnen so nur stellenweise starke 
Anderungen auftreten. Diese Uberlegung gilt fiir jeden Zeitpunkt. 

Der Einfachheit halber wird im weiteren vorausgesetzt, dass sich die Kérner 
bis zu einer beliebig kleinen Grésse zermahlen lassen. Wenn dies nicht zutrifft, 
etwa wegen der Festigkeit der idealen Gitterbereiche, braucht man im folgenden 
nur die untere Grenze des gesamten Korngréssenbereiches von Null auf eine 
minimale Korngrésse zu verdndern und zu beachten, dass e(€, x) nur in dem so 
beschrankten Bereich von Null verschieden sein kann. 


Die Grundgleichung 


Die zeitliche Veranderung 0M(x, t)/0t besteht aus den Ubergangen zweier 
Teilchenarten : 
a) Uberginge der Triimmer der Kérner, die vor der Zerspaltung grésser als 
x + dx waren. Die Masse solcher Triimmer, die pro Zeiteinheit in das be- 
trachtete Intervall (x, x + dx) itbergehen, betragt zur Zeit ¢ 


xm 


‘i e(E, x) M(é, 1) dé. 


1) W. Kiem und A. SmrxKat, Naturwissenschaften 29, 688 (1941). 
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‘b) Ubergange der Triimmer der Kérner, deren eae vor der Zerspaltung im 
Intervall (x, « + dx) lag. Die Masse der Triimmer, die pro Zeiteinheit dieses 
Intervall verlassen, betragt zur Zeit ¢ 


” 
a 
“ 

a 
= 
ie 


a 
ag x 


M(x, i) [ etx, a) da . 


0 


Es gilt daher allgemein 


er + Mx) [ ele, o) da = f o(€, 2) MEE) db. al) 


0 


Diese partielle Integrodifferentialgleichung ist die Grundgleichung der vorlie- 
_ genden Theorie. 

Da beim Mahlen nur kleinere Teilchen aus grésseren entstehen kénnen 
4 («reine Zerkleinerung»), gibt es keine anderen Uberginge als die oben ange- 
_ fiihrten, vorausgesetzt, dass keine Teilchen aus der Miihle verschwinden. Die 
_ Grundgleichung (1) muss deshalb das Prinzip der Erhaltung der gesamten 
Masse erfiillen. Tatsachlich folgt aus (1) 


aie 
) sy | Me) ee 
. 0 


da 


/ dx / e(E, x) M(E, t) dé — / dx i e(x, «) M(x, t) da = 0 


gilt. Da namlich e«(&, «) im Bereiche § < x gleich Null ist, verschwinden die 


Integrale 
i o(E, x) M(é,t) dé, 4 e(E, x) dx 
0 gE 


Deswegen ist 


[asf panes ae fide as t) dé 
a [a Te. “ae x) M(E, t) dx Tala. / e(E, x) M(E, t) dx 


was zu beweisen war. 
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Ermittlung der Funktion ¢(&, x). Experimentelles 


Die Funktion «(&, x) enthalt implizite alle spezifischen Eigenschaften des 
gemahlenen Materials und der Mahlvorrichtung, die mit der Zerkleinerung zu- 
sammenhiingen. Sie ist viel zu kompliziert, um direkt aus Messungen bestimmt 
zu werden. Im folgenden wird deshalb gezeigt, wie sie aus einer fiir Messungen 
geeigneten Annaherung indirekt bestimmt werden kann. 

Der gesamte Korngréssenbereich 0 < x < x,, wird in m Teilbereiche mit 
den Grenzen x) = 0, %4, Xa) +++) %j +++) Xm unterteilt. Dabei wird vorausgesetzt, 
dass 


| PG AG) | < Aa Fie + (2) 


sei. Durch M,,(t) wird die Masse der K6rner bezeichnet, deren Gréssen zur Zeit ¢ 
zwischen x,;_, und x; legen: 


M,(t) = / M(s, t) dx. 


Durch A;, wird weiter die Grésse 


*k xy 
panera eek 

Rea tia 
bezeichnet. Der physikalische Inhalt der Bedingung (2) besteht in folgendem. 
1. Die Teilbereiche sind so klein gegeniiber den entsprechenden Korngrés- 
sen, dass, wenn ein Korn zertriimmert wird, die Masse derjenigen der so ent- 
standenen Triimmer, deren Gréssen im Teilbereich des urspriinglichen Kornes 

bleiben, vernachlassigt werden kann, das heisst 


xy 


i ie / e(E, x) M(E, t) dx = [a / al, x) M(E, t) d& 


%ji-1 *%j-1 


2 ii dx fee x) M(E, 1) dé w0. 


Hj-1 *j-1 


Die ersten zwei Gleichheitszeichen folgen aus der Bedingung der reinen Zer- 
kleinerung 
e(€, x) = 0, wenn f= 2%. 


Hier wurde die Zerkleinerung durch das Reiben und durch Stésse zwischen den 
K6rnern, wodurch von den K6érnern sehr kleine Teilchen abgetrennt werden, 
vernachlassigt. Dies ist zulassig, da das Volumen des Mahlgutes klein gegen- 
iiber dem Volumen der Miihle vorausgesetzt wird und da ausserdem das Mahlen 
mit Mahlkérpern wesentlich schneller vor sich geht. 


- 
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2. Bei der Mittelwertbildung 

© | M&,t) dé = M(*x, 2) (x, — x,_1), 
- 4j-1 

xk 


afecom M(E, t) dé = M(éx' 11 fax fo (Eoavak (ee) Se < uy) 


%h—-1 %j-1 %he—-1  *%i-1 
_ kann wegen der Bedingung (2) 
M (ix, t) » M(‘x’, t) 


- gesetzt werden. 

3 Man integriere nun die Grundgleichung (1) in einem der Teilbereiche, 
4-4 Sx S %x;. Unter Anwendung der Bedingung (2) ergibt sich fiir die ein- 
_ zelnen Glieder 


a i... M(é, t) dé = 5) M(x, t) p fax fot , ae = a Res 14,18): 


k=i+l1 k=1+1 
-1 *k-1 


so dass die Grundgleichung (1) unter der Voraussetzung (2) in ein System von 
_m gewohnlichen Differentialgleichungen 


He + (Stas) Me Y decMy | (webs, et) (3) 


k=i+1 


zerfallt. Die Gesamtheit der Lésungen des Systems (3), die durch die tiblichen 
Integrationsmethoden berechnet werden, stellt schon eine angenaherte Lésung 
der Grundgleichung (1) dar. Die Anfangsbedingungen M;(0) ergeben sich aus 
der bekannten Anfangsverteilung M(x, 0). 

Der zeitliche Verlauf der Gréssen M,(t) kann experimentell leicht verfolgt 
werden, indem das Mahlen in regelmassigen Zeitabstanden unterbrochen, das 
Mahlgut durch Sieben in die Fraktionen M; getrennt und gewogen wird. Ausser- 
dem wird fiir jedes M; (i +1) die reine zeitliche Abnahme verfolgt, indem 
alle Korner grésser als x; jeweils von vornherein abgesiebt werden; daraus er- 
geben sich die spezifischen Zerfallsgeschwindigkeiten in den durch das Absieben 
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vereinfachten Gleichungen (3) 


dM, or ; 
age tales ap > Ary Ab= 2, 3,223) 
k=1 


Aus diesen beiden Messreihen k6énnen alle A; berechnet und die Giiltigkeit des 
Systems (3) gepriift werden. In einer Arbeit von SEDLATSCHEK und Bass’), in 
der die Naherung (3) direkt aus statistischen Uberlegungen abgeleitet und als 
eine selbstandige Theorie dargestellt worden ist, wurde dies durchgefihrt. Es 
ergab sich eine befriedigende Ubereinstimmung der berechneten Resultate mit 
den Messungen. Auf dieser Ubereinstimmung beruht auch die experimentelle 
Bestatigung der vorliegenden Theorie, die ja eine Verfeinerung der ersteren in 
demselben Giiltigkeitsbereich ist, da sie den Grenzfall von (3) fiir unendlich 
' viele, infinitesimale Teilbereiche darstellt. 

Die unter Anwendung der Naherung (3) aus den Messungen ermittelten A; , 
werden nun zur Bestimmung der Funktion e(¢, x) beniitzt, indem man nach 
der Methode der kleinsten Quadrate fordert, dass der Ausdruck 


x4 


y Ss js (%; — %¢_1) — fax f of, x) a + A(x, x) 


t=2 k=1 $ 
X%h-1 %i-1 


ein Minimum annimmt. Fiir praktische Zwecke werden so die Parameter einer 
geeignet gewahlten Interpolationsfunktion bestimmt. Prinzipiell kann «(€, x) 
nach dem Satz von WEIERSTRASS mit beliebiger Genauigkeit durch Polynome 
approximiert werden, wenn man von den mit steigender Anzahl der Teilbe- 
reiche zunehmenden experimentellen Schwierigkeiten absieht. In diesem stren- 
gen Fall verschwindet jeder der m (m — 1)/2 Klammerausdriicke, mit der zu- 
sdtzlichen Bedingung ¢(x, x) = 0; aus dem so entstandenen System linearer 
Gleichungen kénnen [m (m — 1)/2] + 1 Koeffizienten eines Polynoms bestimmt 
werden. 


Lésungen der Grundgleichung 
Die Grundgleichung (1) hat zwei konstante Lésungen. Die triviale Lésung 
bezieht sich auf eine leer laufende Miihle. Fiir den Endzustand gilt?) 
M(x, co) = 8(x) sf M(x, t) dx, (4) 


0 


da eine «reine Zerkleinerung» vor sich geht und da die gesamte Masse konstant 
bleibt. (4) ist tatsichlich eine Lésung, da durch die Substitution in (1) jedes 


1) K. Sepiatscnek und L, Bass, Powder Metal. Bull. 6, 148 (1953). 
2) 6(x) bedeutet die Diracsche Deltafunktion. 
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_ Glied einzeln verschwindet: das erste wegen der Konstanz von (4); das zweite 
_ Glied verschwindet bei x + 0, indem die 6-Funktion verschwindet, bei x = 0 


wegen dex gleichen Integrationsgrenzen. Das Glied auf der rechten Seite von 


_ (1) verschwindet ebenfalls, da 


OS ee) 


fel x) 6() dé =e(0, x) =0 (x <E<x,). 
Bin geschlossener Form lasst sich auch die angenaherte Lésung fiir grosse (das 
_heisst von x,, wenig verschiedene) Korngréssen darstellen. Das Integral auf der 
rechten Seite von (1) kann bei x > x,, wegen dem kleinen Integrationsbereich 


_ vernachlassigt werden, da der Integrand stetig ist. Integration der so verein- 


_ fachten Gleichung (1) ergibt 


=t f ex, a) do 
Metz, t) = M(x, 0)e ° : (5) 


Diese Lésung entspricht der bekannten halbempirischen Formel von Rosin 
und RAMMLER}), die dadurch zum erstenmal theoretisch begriindet und ihr 


_ Giiltigkeitsbereich festgelegt wird. (O. THEIMER?) leitete diese Formel mit den 
_Methoden der statistischen Mechanik ab, allerdings nur fiir den Fall einer 
_zeitunabhangigen Gleichgewichtsverteilung; dabei liess er die Bedingung der 
_«reinen Zerkleinerung» fallen und fiihrte vereinfachte Voraussetzungen tiber 


die physikalische Beschaffenheit der Korner ein.) 
Die Rosin-Rammlersche Formel lautet in der hier beniitzten Ausdrucks- 


- weise 


teense e 6 


wobei ” und F fiir jede Mahldauer neu, empirisch und unabhangig voneinander 
zu bestimmende Konstanten sind, die so implizite die Zeitabhangigkeit der 
Verteilung enthalten. Man kennt aber keine Verkniipfung dieser zwei Kon- 
stanten, obwohl sie offenbar miteinander zusammenhangen miissen; die An- 
fangsverteilung M(x, 0) bleibt ebenfalls unberiicksichtigt, so dass die Formel 
(6) bestenfalls erst nach einer unbestimmten Mahldauer beniitzt werden kann, 
je nachdem, wie ahnlich M(x, 0) einer Verteilung vom Typ (6) ist. 

Die im folgenden angegebene angenaherte Lésung fiir kleine Korngréssen 
sowie die strenge Lisung des Mahlproblems zeigt, dass M(x, ¢) nur im Bereiche 
der grossen Korngréssen in geschlossener Form (mit einem Exponentialterm) 
dargestellt werden kann. Die theoretische Begriindung der Rosin-Rammlerschen 
Formel ist deshalb nur in diesem Bereiche méglich. Die Anwendung von (6) in 
den iibrigen Bereichen wird dadurch erzwungen, dass die Konstanten und F 


1) P. Rosin und E. Ramer, Koll.-Z. 67, 16 (1934). 
2) O. ToEmeER, Koll.-Z. 128, 1 (1952). 
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zu diesem Zwecke verfiigbar gemacht werden, indem ihr Zusammenhang und 
ihre Zeitabhangigkeit unterdriickt werden durch die in jedem Zeitpunkt wieder- 
holte empirische Bestimmung der zwei Konstanten aus dem Verlauf im ge- 
samten Korngrossenbereich. Ausserdem bewirkt diese wiederholte Bestim- 
~ mung das jeweilige Anpassen der Formel (6) an die zeitabhangigen Auswir- 
kungen der Anfangsverteilung, die sonst unberiicksichtigt bleibt. Der Vergleich 
von (5) mit (6) ergibt als entsprechende Ausdriicke 


n(t) x nt) 
M(x,0) <> Ble 


3 nt (7) 
t | e(x, a) de <> Lael 


Die Naherung (5) wird bei grossem ¢ genauer. Wegen der Eigenschaft der 
«reinen Zerkleinerung» wird namlich bei geniigend grossem ¢ 


M(%,,1) = Mite. t) ‘wenn, %.<— 4° 


Je grosser dann ¢ wird, um so starker wird diese Ungleichheit und um so kleiner 
M¢é,t) im Intervall x <& < ~x,,. Mit dieser Verscharfung der Annaherung 
(5) wachst ihr Giiltigkeitsbereich. Daraus erklart sich die bekannte Erscheinung 
bei der Anwendung der Rosin-Rammlerschen Formel, namlich dass sich die 
Anwendbarkeit dieser Formel mit wachsender Mahldauer verbessert. Die For- 
mel (6) nahert sich der Formel (5), soweit es die jeweilige Ahnlichkeit der Aus- 
driicke (7) erlaubt, da die theoretisch unzulassige Anwendung von (6) in immer 
kleineren Bereichen durch die ad hoc getroffene Wahl der zwei Konstanten er- 
zwungen werden muss. 

Bei kleinen Korngréssen kann das zweite Glied auf der linken Seite von (1) 
vernachlassigt werden, da der Integrationsbereich klein wird und der Integrand 
stetig ist. Es bleibt iibrig 

xm 
en = f eles ) MEE, 8 a8. 


x 
Die Integration nach der Zeit ergibt 


xm t 
a 


M(x, 1) — M(x, 0) = [ ag rf e(E, x) M(é, 2) dr. (8) 
x 0 
Die Substitution 
BS teen eR ue oe ote 
transformiert den Integrationsbereich der Variablen €in das Dreieck 


CSS a * Seen 


5 
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; Die vorige Gleichung geht dadurch iiber in 


o 


a t 
M'(x', t) — M'(x',0) = | d&’ | e'(&", x) M(&, 0) de. 
as, 


Der Kern dieser Integralgleichung, e’(é’, x’), ist im Bereiche 
q 


Desist OS PS Sb 


; stetig, wobei b > x,, eine endliche Zahl ist. Dies ist eine Volterrasche Doppel- 
integralgleichung zweiter Art mit der Lésung?) (die Striche werden im folgenden 
-_ weggelassen) 


% t 
M(x, t) = M(x, 0) — | dé ‘e Dy K,(x, t; 1) M(é,0) de, (9) 
6 ae 
wo K;, die entsprechenden iterierten Kerne sind: 
x t 
K,= el, 2), Kix, #61) =—f dy f ely, #) Kialy,2:6,2) de (10) 
‘f t 
| foder allgemeiner, mit grosserer Freiheit in der Berechnung, 


x t 
K(x, t; £0) = fay [ Kyle t; y, 2) Kyjly, 23, 0) dz. (11) 
&é er 


Unter den vorliegenden Bedingungen ist die resultierende unendliche Reihe 
immer konvergent?). 

Da in diesem Falle der Kern, «(¢, x), nicht von der Zeit abhangt, kann jeder 
iterierte Kern K; als Produkt eines rein zeitabhangigen Faktors mit einem 
zeitunabhangigen Faktor Kj dargestellt werden. Die zweite der Gleichungen 
(10) lautet dann 


Klan, t; 6,2) = — Sa f ely, 2) Ki aly, 8) dy, 
§ 


was durch vollstandige Induktion leicht bewiesen werden kann. 
Die strenge Lésung des gesamten Mahlproblems hat die Form 


£3 
—t f ex, a) dx 


M(x, t) =f(x,t)e ° ; (12) 


1) V. VotTERRA, Legons sur les équations intégrales et les équations intégro-différentielles 


(Gauthier-Villars, Paris 1913). 
2) V. VOLTERRA, a. a. O. 


aay 
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Um f(x, 4) zu berechnen, wird (12) in die Grundgleichung (1) eingesetzt. Es 
ergibt sich > 
=t f ex, ot) dor im, -t f eg, ot) dx 
eee = [ee yHene ° dé 


x 


oder, nach Umformung und Integration nach der Zeit, 


E x 
tm [eG a) dx— fe(x, B) ap 


7 pestis serie 0 oo 


Die Gleichung (13) unterscheidet sich von (8) nur um den exponentiellen Faktor 
bei dem Kern. Sie wird deshalb durch dieselbe Transformation wie (8) auf eine 
Volterrasche Gleichung zweiter Art mit einem stetigen Kern, demselben Inte- 
grationsbereich und denselben Konvergenzeigenschaften der Lésung zuriick- 
gefiihrt. Die strenge Ldésung ist daher tatsachlich durch (12) gegeben, wobei 
f(x, t) durch (9) explizite dargestellt ist, wenn dort M(x, ¢) durch f(x, ¢) ersetzt 
wird. Die entsprechenden iterierten Kerne sind durch (10) oder (11) gegeben, 
wenn man eé(é, «) durch 
€ x 
—t le a) da— [etx B) | 
slésx) ex -° : 


ersetzt. 


Summary 


A mathematical theory of milling processes is developed by deriving and 
solving a partial integro-differential equation which describes the time-dependent 
particle size distribution of the milling charge. This basic equation contains a 
characteristic function of both the mill and the material milled, which is conven- 
iently determined from measurements conducted in accordance with an approx- 
imation derived for this purpose from the general theory. This approximation, 
which can as well be used independently is equally convenient for obtaining 
experimental evidence of the validity of the general theory. It has been established 
in a previous paper that this approximation is in agreement with results obtained 
by experiments. 

Several approximate solutions of the basic. equation as well as the rigorous 
one are worked out and an analysis of the semi-empirical Rosin-Rammler formula 
is given from the point of view of the present theory. 


(Eingegangen: 25. September 1953.) 
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Torsionsschwingungen unter Beriicksichtigung 
der Masse und der Dimpfung der elastischen Glieder') 


Von YIAN-NIAN CHEN, Winterthur?) 


In gewissen mechanischen Schwingungssystemen mit langen Wellen, die zu 
_Torsionsschwingungen erregt werden, spielen die Massen und Dampfungen der 
Wellen selbst eine wesentliche Rolle und diirfen daher nicht vernachlassigt 

werden. Dieses Problem stellt sich vor allem bei den Schiffsmaschinenanlagen ; 
sie bestehen aus einer Kolbenmaschine, heute meistens einem Dieselmotor, 
und einem Propeller, die durch eine lange Wellenleitung verbunden sind. In der 
vorliegenden Arbeit werden Verfahren entwickelt, um solche Probleme zu lésen. 
Die Wellenleitung wird im folgenden « Zwischenleitung » genannt, um Verwechs- 
-lungen mit den elastischen Wellen zu vermeiden, die eben in dieser Zwischen- 
leitung hin- und herwandern. 


1. Die Torsionsschwingungen der Zwischenleitung 


In Figur 1 bedeuten: 
@,, R, das Massentragheitsmoment bzw. den Reibungsbeiwert des Motors; 
0,, R, das Massentragheitsmoment bzw. den Reibungsbeiwert des Propellers; 
@,R das Massentragheitsmoment bzw. den Reibungsbeiwert je Langenein- 
heit der Zwischenleitung ; 


is das polare Flachentragheitsmoment des Querschnittes der Zwischen- 
leitung ; 
G den elastischen Gleitmodul der Zwischenleitung ; 


M sinw t bzw. M e'”! das erregende Drehmoment des Motors. 
Der Flachensatz beschreibt die Bewegung pm des Leitungselementes dx 
(Figur 2): 0. OK 


o dx Ot2 = Ox ax aS 


wobei 6 das zu iibertragende Drehmoment und 6p das Dampfungsmoment 
darstellen. Drehmoment 6 und Verdrehung @ sind durch die Beziehung aus 
der Elastizitatslehre verkniipft : 
0 
NM — 4 G ara . 


1) Die Anregung zu der vorliegenden Arbeit verdanke ich Herrn A. KLEINER, Chefingenieur 
der Forschungsabteilung fiir Dynamik, Gebriider Sulzer AG., Winterthur. 
2) Gebriider Sulzer AG. 
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Das Dampfungsmoment J¢p wirkt immer der Bewegung entgegen, also in 
Gegenphase mit der Geschwindigkeit, daher 


op 


Damit wird | 
he RAPE Mea 5 BG AR Ry Beas (1) 
ot? Oo ot o Ox? 
oe 
OM 
on 


SULZER 
74527 


SULZER 
74526 


Figur 1 Figur 2 
Zweimassensystem mit elastischer Leitungselement. 
Zwischenleitung. 


Die partielle Differentialgleichung wird nach D’ALEMBERT durch die Funktion 
befriedigt : 
p=ff{e(ct+ x)}+ gfe (ct—x)} [e= Konstante], 


welche zwei Wanderwellen darstellt, die sich in entgegengesetzten Richtungen 
mit der Geschwindigkeit c durch die Leitung fortpflanzen. Einsetzen dieses 
Ausdruckes 9 in die partielle Differentialgleichung ergibt 


e (2 2) tfe(ert x} + e"{elet—x)}] 
a F clf'{e (ct + x) }+ g'{e(ct — 2) t1=0., 


wenn aos als Abkiirzung von att... di san f usw. eingefiihrt wird. Diese 
Gleichung kann noch weiter vereinfacht werden, falls folgende Beziehung 
besteht 


f'{e (ct+x)}+ e"fe (ct — a) = d[f'{e (ct + x) }+ g'{e (ct — x) ¥] 


mit 6 als Konstante. Man erhialt dann 


[r"{e(ct+a)}+ etfelct— a} {octet +5 Acoli 
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oder 


gies 


éc% — + bh =, 0\, (2) 
woraus sich die noch unbekannte Geschwindigkeit c ermitteln lasst. 

Diese Vereinfachung wird offenbar durch eine Exponentialfunktion fiir Y 
-ermoglicht, die dann die Lésung bildet, namlich 


Rid i (ct+x+ a’) w/c t(ct—%*+ B’) w/c 
= Ave + Be 
oder 
Q= vA eg? (@t+o x/c+a) be B et @t-o aIc+ B) 


worin @ die Kreisfrequenz des erregenden Drehmomentes und A, B, «, f die 
Integrationskonstanten sind. Es wird spater noch gezeigt, dass die Denne die 
aufzustellenden Randbedingungen befriedigt. Die Konstanten ¢ und 6 in der 
'Gleichung (2) sind dann 


c=0 


Ra, gO = 
C é o 
oder 


= eam = 


Der Ausdruck R/(@ #) ist immer sehr klein gegen 1, darum lasst sich die Formel 
der Geschwindigkeit vereinfachen zu 


pelt Rk? RAO Nate 
owe, (1 ae 2wu% a sorgr) (1+ sir} 


Came te a3 
Die Fortpflanzungsgeschwindigkeit c erscheint in komplexer Form infolge der 
. Wirkung der Dampfung R. Der reziproke Wert von c: 


mit 


+ die Gleichung (2) wird 
i 
: 
, 


1 : 
— = es (a ms of b) 
-enthalt zwei Gréssen 
kR? R 
a= (1+4—793) und See 
-wovon a nur wenig von 1 verschieden und 0 viel kleiner als 1 und hauptsachlich 
von der Dampfung R/(2 w #) bestimmt ist. 


| 
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Die Lésung fiir y wird: 
p(x t) = A eth ealer ei (wttawxler +a) as B ge bw aler et (wt-awx|cr+ B) ; 
Wenn der Ausdruck Impedanz z: 
: ba 1a 
Z=Uu+tvV= 3 Pigs 
eingefiihrt wird, wird der Winkelausschlag y am Ort x zur Zeit ¢: 
p(x, t) ee et ei(wituxrta) -} Be-%* eg? (a t-vx+ B) : (3) 


Wegen der Dampfung, die sich im reellen Teil « der Impedanz z ausdriickt, ist 
es nicht méglich, den Winkelausschlag auf eine reine Sinusfunktion zuriick- 
zufiihren. 


2. Das Zweimassensystem mit elastischer Zwischenleitung 


An beiden Enden der elastischen Zwischenleitung stehen je eine Masse mit 
Tragheits- und Reibungsmoment (Figur 1). Auf der einen Masse 9, wirkt noch 
eine periodische Erregung M e‘®*. 

2.1. Randbedingung an der Stelle x = 1 (Motorseite) 

Die Randbedingung lautet 


0% 0 0 ne 
0, a + Rt eG (52), = Me Ay 


Nach Einsetzen der Gleichung (3) ergibt sich 


[-0* 0, +40 KR + Jp G = (b af ia)h A eb @ller pi (wttae ler +a) 
ae {-0* 0, +740 R, = Jy G = (b +4 a)| B e7 eter gilwt-aeller+p) __ jy pivt : 


Werden die Impedanz der Masse: Z, = —p, +i 9,= —w?0,+i@ Ry und 
der Profilwert der Zwischenleitung: k = J, G eingefiihrt, so erhalt man nach 
Aufheben der Glieder e*”* auf den beiden Seiten der Gleichung 


Ae! (Z, + kz) e+ Be (Z,—k2e = M. (4) 


2.2. Randbedingung an der Stelle x = 0 (Propellerseite) 
Die Bedingung ist hier 


0° Op 0 
Pee ec Koa alae (52), = 0 
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Einsetzen der Gleichung (3) ergibt 

| A ef (Z,— hz) + Be? (Z,+k2z) =0 (5) 
mit 
; Zo = —-07 Oo +10 Ry =—Pyo tig. 


_ Aus den beiden simultanen Gleichungen (4) und (5) folgen die beiden Unbe- 
kannten A e** und Be‘? zu 


Ac*™— M(Z.+k2){(Z, + kz) (Zot 2) ce! — (Z,—h2) (Zy— kz) } 1, 

‘Bei? = -M(Z,—k 2) {(Z1 + h2) (Zot kz) — (Z,— kz) (Z)— kale #}}, 

‘und der Winkelausschlag : 

| (x, t) = M{ (Z, + h 2) e¥* eMltt?#) — (Z, — hz) e-** efoto) A 
x {(Z, + kz) (Zo +h) e*— (Z,— kz) (Zy— ka) ec}, 


-  Somit ist erwiesen, dass der Lésungsansatz der Exponentialfunktion die 
-Randbedingungen. erfiillt. 

__ Die ganze abgeleitete Beziehung lasst sich im folgenden Abschnitt mit Vek- 
toren noch viel deutlicher darstellen. 


3. Vektorielle Darstellungen 


Die Randbedingung auf der Propellerseite, Gleichung (5), lasst sich unter 
Einfiihrung der Vektoren A= A e’* und B = Be’? umformen zu!) 


Z) (A+ B) —k2(A—8)=0. (7) 
Zar vektoriellen Darstellung dieser Gleichung nimmt man zuniachst einen Ein- 
heitsvektor Poet. Ak 
arate pA Bh. 


an und bildet aus ihm den Vektor Z, (A) + By) (vgl. Figur 3). Mit diesem als 
Durchmesser zeichnet man einen Halbkreis, dann ist die vom Anfangspunkt O 
unter Winkel « [ctg « = b/a = R/(2 w #)] gezogene Sehne der Vektor ku (Ay — By), 
woraus sich der Vektor (Ay — By) bestimmen lasst. Aus (Ay + By) und (Ay — By) 
k6nnen die beiden Vektoren A, und By konstruiert werden (siehe Figur 3). 

Die Randbedingung auf der Motorseite aus Gleichung (4) 


{Z,+kz}Ac!+{Z,—kz}Be"=M (8) 
kann dann von den beiden Vektoren A, und By ausgehend dargestellt werden. 


1) O, Lurz, Resonanzschwingungen in den Rohrleitungen von Kolbenmaschinen, Ber. Labor. 
Verbrennungskraftmaschinen, Technische Hochschule Stuttgart, Heft 3 (1934). 
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79, (Ay+B,) 
SULZER -f, (A,*8,) 0 A,B =1 
74528 

Figur 3 


Bestimmung der Vektoren A) und By aus der Randbedingung auf der Propellerseite. 


SULZER 
74529 


Figur 4 
Bestimmung des Vektors My aus den beiden Vektoren Ap und By. 


ZAMP 
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: In Figur 4 wird zunachst aus Ay der Vektor (Z, + k 2) Ap konstruiert und dann 
-skalar mit e“! multipliziert. Nach Drehen des neuen Vektors um den Winkel v / 
kommt der Vektor (Z, + k 2) Ay e*! zustande. Analog bildet man aus Vektor By 

der Reihenfolge mach die Vektoren (Z,—k2z) By, (2,—k 2) Bo enh, 


_ (Z,—k2) Bye-**. Aus den beiden Vektoren (Z, + #2) Age*' und (Z,—k2) Bye-*! 


i ae Ux 


Figur 5 8 
Bestimmung des Winkelausschlages fiir eine beliebige Stelle zur Zeit ¢ = 0. 


 lasst sich die Resultierende M, bilden. Nach Gleichung (8) besteht dann die 


Beziehung im B 

Ni ee ve Fae ee 
oder umgekehrt 

M M 

A as Ay M,’ B = Bo M, . (9) 


Somit sind die beiden Vektoren A und B bekannt. 

Fiir irgendeine Stelle x lasst sich der Winkelausschlag y zur Zeit ¢ gemass 
- Gleichung (3) wie folgt ermitteln. In Figur 5 geht man von den Vektoren A 
und B aus und bildet Ae“” und Be-“®, diese werden um den Winkel ++ v x bzw. 
—vx gedreht und deren Resultierende gebildet. Dann ist der neue Vektor 


gleich dem Winkelausschlag an der Stelle x zur Zeit ¢ = 0: 
p(x,f=0) = Ae + Bee. 


Nach weiterer Drehung dieses Vektors um den Winkel w¢ erhalt man den 
Ausschlag zur Zeit t: p(x, t) = p(x, 0) e*?! (vgl. Figur 6). 

Da das erregende Drehmoment M sinw ¢ tatsaichlich nur dem imagindren 
Teil von M e‘” entspricht, ist der wirkliche Ausschlag auch nur’ durch den 
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imaginaren Teil von g(x, t) dargestellt, also nur durch dessen vertikale Kom- : 
ponente, wie in Figur 6 gezeigt. Der Winkelausschlag an der Stelle x - 0 zur 
Zeit t: (0, t) == y(0) e*?* baw. p(0) sinw ¢ lasst sich auch zeichnen. Der Winkel y 


zwischen den Vektoren g(x, 0) und g(0) stellt ihre Phasenverschiebung dar. Im 
dampfungsfreien Fall ist y = 0 oder z. 

Fallen die beiden Vektoren (Z, + kz) Ay e?! und (Z, — kz) Bye~*' in Figur 4 
in gleiche Richtung, ist My am grossten, dann sind A, B und damit g am klein- 
sten. Das heisst, das System wird sehr schwach zum Schwingen erregt. 


~(S 


SULZER SULZER 
7453) 74532 


Figur 6 Figur 7 
Die Winkelausschlage q(x, t) und (0, 2). Die Restgliedkurve Al/l = f(a). 


Ein Resonanzfall liegt aber vor, wenn die beiden erwahnten Vektorkompo- 
nenten von M, entgegengesetzt liegen. Somit lasst sich die Eigenfrequenz auf 
folgendem Weg bestimmen. Zunachst wird ein passendes w als Kreisfrequenz 
gewahlt und diese beiden Vektoren konstruiert. Man ermittelt den nétigen 
Winkel v Al, um sie in entgegengesetzte Richtungen zu bringen. Fiir eine 
gewisse Lange / + Al der Leitung des Schwingungssystemes ist die gewahlte 
Kreisfrequenz w bereits die Eigenkreisfrequenz. Fiir eine Lange / wird die 
Eigenkreisfrequenz 

‘1+ Al\n 
WO, = @ ( —— ) 


l 


sein, wobei der Exponent » zwischen 1/2 und 1 liegt und noch abzuschitzen ist. 
Sind die Massen 0, und @, iitberwiegend gross, so ist » ungefahr gleich 1/2, 
sonst wird bei sehr langer Zwischenleitung nahe 1 liegen. Nach Abschatzung 
von # erhalt man eine neue Kreisfrequenz und wiederholt mit ihr den ganzen 
Vorgang nochmals, und so weiter. Die so gefundene Grésse Al/l wird tiber w 
zu einer Kurve aufgetragen, deren Schnittpunkt mit der w-Achse die Eigen- 
frequenz darstellt (siehe Figur 7). 


3 V, 1954 Torsionsschwingungen 301 


| 4. Mehrmassensystem mit Erregung an nur einer Endmasse 


Figur 8 stellt ein (n + 1)-Massensystem dar. Ein Erregungsmoment M sin wt 
wirkt auf die Endmasse x. Fiir eine beliebige Masse 7 mit Massentragheits- 
moment @, gilt pes cect 

0 0 
oe, 2 -Ot2- LN og pe Ot eG =e! — Jon G PE = 0 
mit 


= 4} EtIe gOt aerated) B; 7 H5* gt(wt—vjx+ Bj) 


fiir den linken Leitungsteil 7; 


4 
; 
@, 9, @, a 2: On Msin wt 
i 
: 
_ 
: 
: 


re 


SULZER 
74533 


yj : 
Figur 8 


(nm + 1)-Massensystem mit Erregung an der Endmasse 1. 


yee A, etk* ¢ t(@ t+ V7_ x + ae) A. B, md eiot—vex+ Br) 
fiir den rechten Leitungsteil ;; 


R, = Reibungsbeiwert der Masse 7; 
_ Jp; = polares Querschnittstragheitsmoment im Leitungsteil 7; 

pk = polares Querschnittstragheitsmoment im Leitungsteil k. 
_ Aus gy; und g; werden die Ableitungen nach x und ¢ gebildet und in die obige 
_ Differentialgleichung eingesetzt. Unter Beriicksichtigung der Kontinuitatsbe- 


dingung y = 9; = 9; bei x =/1, erhalt man 


1 "2 aoe ate ia 
oh, te, |? Z{ (Z; + k z;) A; erst +- (Z; ad k 2;) B; é 141} = An é kl = B, é al j 
“und 
. A, eral 4B, ett — A, erkls 4 B, ge thy ; 


wobei z, =“, —1v;, die konjugierte Impedanz von 2, ist. Die beiden Glei- 
- chungen stellen die Randbedingung der Masse 7 dar. 

Waren A, e*/47 und B,e~*'3 bekannt, dann kénnten A, e**) und B,e~**4 aus 
dieser Pondbetngine jy seagaat werden. In Figur 9 wird gezeigt, dass aus 
A, e*i4 und B;e~*#'i die Vektoren (Z; + ; 2;) A; e747 und (Z; — hk; 2;) B; e779 

Beh bilden ficsehe Multipliziert man deren F Aaleonde (V) mit 


1 1 
SU DLW St - a De Kae 
k ky, | 2% |? 
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so erhalt man zwei neue Vektoren und bildet deren Summe A,,e**/4 — B, e*k'i, Da 
anderseits aus der Kontinuitatsgleichung die Resultierende der Vektoren A, e*kly 
und B, e~**! bekannt ist, kénnen diese beiden nach der Darstellung in Figur 9 
ermittelt werden. 

Physikalisch bedeutet das Ergebnis der mathematischen Entwicklung, dass 
die beiden Wellen B, und A, im linken Leitungsteil 7 tiber die Masse 7 in den bzw. 
aus dem rechten Teil # wandern, unter gleichzeitiger Anderung ihrer Amplitu- 
den auf B, bzw. A,. 


SULZER 
wae rata) 74534 
7 


Figur 9 


Konstruktion der Vektoren A,, e*k!li und Bye “#k!5 aus jenen A; e*3!7 und Boe *3'3, Nach den Buch- 
staben A und B in der Figur ist jeweils der Exponentialausdruck e** der Einfachheit halber 
weggelassen. 


Fir die Endmasse 0 ist die Randbedingung nach Gleichung (7): 
Zo (Ay + By) — ky % (Ar — By) = 0. 


Die Berechnung kann somit von dieser Masse aus anfangen. Man nimmt den 
Einheitsvektor oan 
(Aro + Bio) = | a 
an und ermittelt die Vektoren Aj) und By, nach Figur 3. Bis diese beiden Wellen 
die Masse 1 erreichen, sind sie zu Ay) e*" und By) e~*" geworden und kénnen 
wie in Figur 5 konstruiert werden. Uber die Masse 1 verwandeln sie sich zu 
Ago e* und Boy e~*4, deren Ermittlung bereits in Figur 9 wiedergegeben ist. 
Sie wandern weiter fort tiber die Massen 2, 3, 4, ... bis zur letzten », an welcher 
die Beziehung der Gleichung (8) gilt. Durch eine graphische Berechnung wie 
in Figur 4 lasst sich der Vektor My bestimmen. Nach Gleichung (9) sind dann 

die Wellen in einem beliebigen Leitungsteil 7 bekannt, namlich 


M 


a ie io MM,’ 


M 


_ Vol. V, 1954 Torsionsschwingungen 303 , 


y 

em der Winkelausschlag g an dieser Stelle nach Gleichung (3) bestimmt 
werden kann. 

: 

5. Dreimassensystem mit Erregung an beiden Endmassen 


Wenn in einem Mehrmassensystem mehr als eine Masse durch dussere Mo- 
“mente erregt wird, dann kann man keine Lésung in analytischer Form erbrin- 

gen. Es miissen die zeitlichen Fortpflanzungswege der von den erregenden 
_ Impulsen ausgelésten Wellen verfolgt werden, bis stationdrer Zustand erreicht 
a wird, das heisst, bis sich tiberall stehende Wellen bilden. Diese Fortpflanzungs- 
-wege im (x, t)-Diagramm heissen Charakteristiken. 


aa 2 


> 


ee ee ee ee eee 


Jisin ot wl sin(wk-d) 
i / 2 i 
© Pi Osun 
eee 
Figur 10 


Dreimassensystem mit Erregung an beiden Endmassen. 


In Figur 10 wird ein Dreimassensystem gezeigt, an dessen beiden End- 
massen die erregenden Drehmomente wirken. Man fangt mit dem Anlaufvor- 
_ gang an. Von den beiden Endmassen 0 und 2 gehen die erregenden Impulse als 


~ wandernde Wellen nach der Masse 1 hin. Die von der Masse 0 ausgehende Welle 
y : : 

_ besteht jetzt aus dem hingehenden Teil allein, also 

. p= a et (mtu *+ Bi) 

und wird durch die Randbedingung: 

i 

i 0°y op Gay pte 

Oo a2 + Boj — Jor G oy = Me’ 


"bestimmt. Nach Einsetzen des Ausdruckes ¢ erhalt man 
By = —M (—Z, — ky 2,)7' . (10) 


- Diese Beziehung stellt zugleich den momentanen maximalen Schwingungsaus- 
schlag der Masse 0, erregt durch das Moment M e'®’, dar. Es sei angenommen, 
_ dir Welle B, erreiche die Masse 1 frither als die von der Masse 2 ausgehende und 
_ besitze bei ihrer Ankunft in 1 den Wert: 


p= B, en ah gt(ot-% 1+ Bi) . 


] 


Diese Welle kann nur teilweise iiber die Masse 1 auf den Leitungsteil 2 iiber- 
gehen, wahrend der tibrige Teil in den Leitungsteil 1 zuriickreflektiert wird. 
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Diese beiden Teile sind 


B, en taht et (wt veh + Bs) und As eth eilotrih +a) 


und miissen mit der urspriinglichen Welle folgender Randbedingung an der 
Masse 1 geniigen: 


OP2 
6, ° se 7 iio P+ Jo1 G = Sox G ae eee 0. 
Man setzt 
%1= ae etal eilmtt ul + a) ae rm eo tah et (etl + Br) 
i= 
und 
Gs, = Bee” sl, pi(wt— V2 1, + Bs) 


in die Differentialgleichung ein und erhalt 
(Z, + &, 2) Ape*™ + (Z, — hy m) Be 8 S Rye, (—By en th) 
Wird das Glied B, e~# aus dieser Gleichung und der Kontinuitatsgleichung 
A, e%h + Beh = Beh 


eliminiert, so ergibt sich 


A pile, See —hk, 21 + Re 2g 
? = 


ey4; 
= pees apy eR ay ; (11) 


Figur 11 veranschaulicht die Konstruktion der beiden Wellen A, e*" und Be" 

aus dem erregenden Drehmoment M. Die Nummern im Kreis verfolgen den 

Vorgang. 

1. Gemass Gleichung (10) ergibt sich B, aus —M. 

2. Aus B, wird B, e~*" gezeichnet. 

3. Aus B, e~* geht (Z, — hy 2, + Rp 2) B, e~*" hervor. 

4. Gemass Gleichung (11) ist der Vektor (—Z, — hk, 2, — ky 25) A, e* bekannt. 
Damit k6nnen seine beiden Komponenten 


(py — hy Uy — hy Me) Aye®"™ und 4 (—g, — A, 0, — Ry vQ) A, e*** 
konstruiert werden, woraus A, e*” bestimmbar ist. 


5. Die Vektoren A, e*” und B, e~*" addieren sich schliesslich zu B, e~*". 


Die Welle B, e ** wandert nach der Masse 2 fort, trifft unterwegs die von 
dieser ausgehende Welle 


Meee ee ae - ; 
V> = Se SIe gmat potty 2) [vgl. Gleichung (10)], 
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kreuzt sie, erreicht dann die Masse 2um die Zeit t und wird als Welle A, er! 
reflektiert. Die Randbedingung ist hierbei 


02 
0, rs aR 


0 
2 - t Jon @ =e = M efter 
mit 


twt —2 7 
p= A,e** 6° + Boe ** ef?! x=], und t=. 


Figur 11 


Konstruktion der beiden Vektoren A, e” 4 und Bea h 


| 


(—Z,+ ko 5} B, e7 #3 aN: Me-*® ie (Zi. ae ky 29) A, ae : 


*l aus M e~*? 


Der Vorgang der Konstruktion von A, é 
Figur 12 gezeigt. 


1. Aus B, e~*! (von Figur 11 iibernommen) wird B, e~*" konstruiert. 


ee ee ee ee eee 


. _* => 


wird. 


— 


_ ZAMP V/20 


Nach Einfiihrung des Ausdruckes @ in die Differentialgleichung ergibt sich 


3. Aus diesem Vektor findet man (—p, + hy uy) Ay e**" und weiter A, e#” 
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eier 


aus dem erregenden Drehmoment M. 


(12) 


und B,e *" wird in 


2. Gleichung (12) lasst sich dann darstellen, womit (Z, + hy 25) A, e*'= bekannt 
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Aye Zl: 


SULZER 
74537 


Figur 12 


: ‘ole -4 — aly 
Kxonstruktion des Vektors A, eta! aus Me a Td Boe. 


Figur 13a stellt die Fortpflanzungswege oder Charakteristiken der Wellen 
im (x, t)-Diagramm dar. Die Wege sind in jedem Leitungsteil gerade Linien, 
ihre Neigung bedeutet die Fortpflanzungsgeschwindigkeit, namlich 


tea; = (y mit ge l2. 


Die Wege erfahren eine Knickung, wo sie durch die Mittelmasse 1 gehen. In 
Figur 136 werden die Amplituden der Wellen in logarithmischem MaBstab 
dargestellt, zum Beispiel fiir Welle B: 


In (Be-“*) = In B+ Ine“**=n B—4z., 


Diese Kurven sind wieder Gerade mit Neigungen tg 8 = — wu fir die hinlaufen- 
den und tg 6 = +w fiir die riicklaufenden Wellen. Fiir irgendeine Stelle x 
kann somit die Wellenamplitude abgelesen werden. Sie nimmt mit zunehmen- 
dem Wege infolge der Dampfung b ab. 

Die oben behandelten Wellen B,, A,, By und A, sind bereits in Figur 13 auf- 
gezeichnet. 

Nun kommt man auf die von der Masse 2 ausgeléste Welle Aj zuriick. Beim 
Erreichen der Masse 1 hat sie den Wert 


aM at gah giat 

—Zy— hy ds ; 
Sie kreuzt hier die von der Masse 0 kommende Welle B, e~*" e*®* und wandert 
auf den Leitungsteil 1 itber, mit veranderter Starke und unter gleichzeitiger 
Reflexion. Die beiden neu entstehenden Wellen Aj und By lassen sich aus der 
Beziehung : 


= Ai enh eit 4. B, ew tats eivt 
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Ee ee Re rl ee 
: . iy 


i ati 


Ce ee Ce ENT ee EE eee yee ren ee 


SULZER 
74538 


Figur 13 


a Wellenplan im (x, ¢)-Diagramm. } Abhangigkeit der Amplituden der Wellen mit dem Weg x. 


und _m 


ae! way siee S et giv t-9) a3 B. eo ahs eivt 
2 


| gemiss der Randbedingung an der Masse 1 aus den Wellen B, und A, er- 


; mitteln. 


Auf diese Weise kénnen die weiteren Wellen verfolgt werden, bis tiberall 


-stehende Wellen entstehen. 


6. Schwingungssystem mit Verzweigung 


Figur 14 zeigt ein Dreimassensystem mit einer angehangten Masse P an der 


Mittelmasse O durch eine Leitung 3. Man stelle sich etwa vor, die beiden End- 
-massen M und M’ seien zwei Dieselmotoren, welche durch ein gemeinsames 
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Getriebe O den Propeller P antreiben. An der Masse O miissen folgende Rand- 
ai: oe werden: 


0 Oy e. 
oe, © om “i ee aa I ATE: \ — JxoG <i Ive G aoe a=i() bei x = L, 
und 
V1 = Po = Ps 
mit 


G3 = A eats eivt ae B, eh eit (j thes i 2 3) : 


(2) ®, Oy 


(a) SULZER 
D 74539 
Figur 14 


Dreimassensystem mit einer angehangten Masse. 


Daraus ergibt sich 
(Zo + Ry %) Ape” + (Zy — hy %) Bre" 
= hy %y (A, 78 — By eh) + Ry 2q (Ay 62 — By e 8") , 


Aus den drei nach der Masse O wandernden Wellen B,, A, und A, kénnen die 
reflektierten Wellen A,, B, und B, gemass dieser Gleichung bestimmt werden. 
Man fangt ebenfalls mit dem Anlaufvorgang an. Es gehen jetzt die Wellen B, 
und A, von der Masse M bzw. M’ nach derjenigen O hin, wahrend die Welle A, 
in der Leitung 3 noch fehlt. Erst nachdem eine der beiden vorhergenannten 
Wellen die Masse O erreicht hat, entstehen im Leitungsteil 1,2 und 3 die 
Wellen A,, B, bzw. B;. Wenn die Welle B, am Propeller P ankommt, bildet 
sich als ihre Reflexion die Welle A;. Die weitere Behandlung der Wellen erfolgt 
sinngemass nach demselben Prinzip wie im vorigen Abschnitt. 


7. Zerlegung eines Schwingungssystems in einzelne Teilsysteme 


Eine Schiffsdieselanlage wird hier als Beispiel gewahlt. Bei der Berechnung 
der Torsionsschwingungen der Kurbelwelle ist es iiblich, ihre Massen und Ela- 
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tizitaten getrennt zu behandeln. Die Kurbelwelle wird namlich in einzelne 
individuelle Massen zerlegt und die dazwischenliegenden elastischen Teile als 
massenlos und dampfungsfrei angesehen. Anders ist es in einer langen Leitung, 
"welche die Dieselmaschine mit dem Propeller verbindet. Ihre Masse muss dann 
“als kontinuierlich verteilt und mit Elastizitat und Dampfung behaftet behandelt 
werden. Fiir sie kommt die in der vorliegenden Abhandlung entwickelte Methode 
_zar Anwendung und fir die Kurbelwelle die wesentlich einfacheren wohlbe- 
5 "kannten Verfahren von HoizeER, BRABEC und andern. Es ist aber schwer, in einem 
" einzigen System zwei Teile nach zwei verschiedenen Methoden zu behandeln. 
' Man sieht sich daher vor die Aufgabe gestellt, die beiden Teile zu trennen. 
_ In Figur 15 stellen 0, bis @,, die Dieselmaschine dar, 9,,,, das Schwungrad 
ind @,, mm den Propeller. Zwischen den beiden letzteren liegt die lange Leitung 
mit lauter angebauten Massen @,, , 9, On +3, «++» On + m—1. Das Schwungrad trennt 
das System also in zwei Teile verschiedener Natur. Fiir die 1 + m + 1 Massen 
_ lassen sich folgende Bewegungsgleichungen anschreiben : 


; 
% (Se), + Bo (Sr), — bs eral et (0) 
, (32), Pes (3°), + Jn © (324, Talon Gear Me, (1) 
J 

f I Oe er eee re, oie cate G8 lone eho ads, Nd wu aito (wae Spe see) Fis we algh metas S200 
Be, ) 

: NYE em Oban (n) 
. re 42) fons G(s), = een, 

: 

Ons (3), =e Recey (3)... | (n ca 1) 
BiG) Jone (24), = 0, 

| O,, +2 (30), + Rie (3 lee: | (n +2) 
) 

ot (Sets) SAG (“Set4) 0, | 
ED ee 
ee tr. (7)... ae 
can — 
A RES ea ee Salant nal (“") sig 


0*y oy OP n+ A ng; 
Onsm i ea é Rn +m cle +m x Jonem G Ox In +m , b m te) 
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wobei M, e’(!+%) das erregende Moment auf die Masse 0; (j= 1, 2,3, ..., ) : 
bedeutet und y seine Ordnungszahl. Spaltet man die Gleichung (7 + 1) in die 
zwei folgenden ] 


tas (Gr) + Rha (GP), tdowaae (“Set), =0, (+a) 
0 


0 
rt (32) 4+ Rat, (GP), — Jonze (4), = 0, @ ei 


wobei 
— * . i 
Ot ., + OF4, =0,4, und RE, + Riti= faa 


sind, so lasst sich das ganze Gleichungssystem in zwei Gruppen zerlegen. Die 
erste Gruppe umfasst die Gleichungen (0), (1), ..., (w) und (m + 1)*, die zweite 
die Gleichungen (7 + 1)**, (w+ 2),..., (7+ m). Es ist ersichtlich, dass die 
erste Gruppe eine einfache Dieselmaschine mit dem Schwungrad O%,, dar- 
stellt und die zweite das m-Massen-System O**,, 9,15, ---, On 4 m- In diesem 
Sinne kann man sich das urspriingliche System in die zwei entsprechenden 
Teilsysteme (siehe Figur 16) zerlegt denken. Damit die Schwingungen dann 
noch unverandert bleiben, muss die Aufteilung von 9, ,, und R,, ,, nach einer 
eindeutigen Regel erfolgen. 

Das erste Teilsystem fiihrt eine erzwungene Schwingung aus, welche durch 
M, e'**°s) erregt wird. Somit ist w die Kreisfrequenz dieser erzwungenen 


n+l N+2 n+s n+n-t nem 


SULZER 
74540 


Figur 16 


Figur 15. — Das (n + 1+ m)-Massen-System stellt eine Schiffsdieselanlage dar; @y bis Q,, die Diesel- 
maschine, @,,,, das Schwungrad, Ose den Propeller und die Leitung zwischen @ 
und @, |, die Zwischenleitung. 

Figur 16.— Das System wird in zwei Gruppen zetlegt: Oy bis ON die Dieselmaschine mit einem 
Schwungrad, und oer, bis Oran die Zwischenleitung mit m Massen. 


n+1 
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Die Baifoabe besteht also darin, eine Masse O**, und ihre Dampfung R**, zu 
bestimmen, mit welcher das zweite System eine Eigenkreisfrequenz w erwerben 
kann. Ist dies getan, so kann das Teilschwungrad 0% ,, fiir das erste System 
a O77 0,4, O77, und dessen Dampfung aus R*,, = R,,,— Rit; 
stimmt werden. Da dieses Teilsystem sich mit irgendeinem der iiblichen Ver- 


fahren behandeln lasst, ist der Schwingungsausschlag von @* , , leicht zu be- 


—- 
< SULZER 
74543 


Figur 17 
Die Zwischenleitung. Von ihrem rechten Ende aus wird x gerechnet. 


f 
sstimmen. Darauf wird nicht weiter eingegangen. Dieser Ausschlag wird auf die 
“Masse O**, des zweiten Systemes iibertragen, woinit die Randbedingung fiir 
dieses bestimmt ist. 

Ss Im folgenden wird die Bestimmung von O**, und R**, erdrtert. Um die 
‘Berechnung zu erleichtern, fangt man bei der Endmasse 0,,4m an, von wo 
aus langs x gerechnet wird (siehe Figur 17), genau so wie im Abschnitt 4 be- 
‘schrieben. Nach Annahme des Einheitsvektors Aj: m), + Bi+m), an der Masse 
nm k6nnen iiber die Massen O,,, m—1 On+m—g «++» On4g bis zur Masse O74, 
die Vektoren Aj, 2), emm+2'n+s und Biya), 6 Mt? n+: bestimmt werden. Fiir die 


“Masse O**, gilt gemass Gleichung (7) die Randbedingung 
oh einseln+s as B., 12), -Sntaln4 ') 


z; 1 -~Zn+el 
=. Zn+2'n =< n+2o'n+i 
<a Ry io 2a x2 (Ain 2), & i" of B, 2), ? : ) 


; 
4 
} 
B 
J 
4 
, 
j 


‘oder, abgekiirzt geschrieben, 


| —Z**, (A+B) = Ry. 2412 (A— 8)- 


Figur 18 stellt diese Gleichung dar. Man bildet aus A und B die Vektoren A+ B 
und A— B. Aus dem letzteren lasst sich der Vektor ky +9 2n 42 (A— 8) konstru- 
jeren, von dessen Spitze eine Senkrechte auf den Vektor A+ B gerichtet wird. 


7 
- 
j 
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Dann sind die beiden Strecken 


9 = PAE, (A+B) = 07 OFF, (A+ 8) 
und 
i = —i g**, (A+B) = —i w R**, (A+ B) 
und damit 


t= F{ot(A+B)}* und Ret, = sW{o (A+ B)}. 


n+1 n+1 


Die Aufgabe der Zerlegung der Masse @,,., und ihrer Dampfung R,,,, ist 
damit gelost. 


SULZER 


ee 7454? 


Figur 18 


n+1 
se—= 5) 


Fs Ke : - Ingol 
Bestimmung der Gréssen Ox und R**_ aus den beiden Vektoren A = Ain 2), & he 2 eee 


ge inteln+s 
(n+ 2), ‘ 


@*x*, kann negativ sein, wenn y und A+B entgegengerichtet sind. Da- 
gegen bleibt die Dampfung R**, immer negativ. Der Ausdruck 7 w R**, @nay 
stellt dann das vom ersten auf das zweite System iibertragene Drehmoment 
dar. Dieses dient dazu, um die Energiebilanz im zweiten System aufrechtzu- 
erhalten, da seine Dampfungen mechanische Energie fortwahrend vernichten. 

Das Ergebnis kann noch anders gedeutet werden: Der Teil rechts vom 
Schwungrad lasst sich vollstandig auf eine einzige Masse mit einer einzigen 
Dampfung reduzieren. 


8. Torsionsschwingungen ohne Dampfung 


Eine Nockenwelle tragt sehr viel kleine Schwungmassen. Um ihre Eigen- 
frequenz zu bestimmen, geniigt es in den meisten Fallen, die Dampfungen 
ausser acht zu lassen. In Figur 19 ist eine Nockenwelle mit Antrieb in der Mitte 
dargestellt. Unter Weglassen ihrer Dampfungsglieder sind die Randbedingungen 
fiir x = 0 gemass Gleichung (7): 


Jor G — (Ay — B,) = 0 
1 
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nd fiir eine beliebige Masse 7 nach Abschnitt 4: 


sees |(; Sats ) Ave laleg . w O;C; in 

- 1 af ee” Kt haa (- ae SAID ) Bye iw@ljic \ 

4 re ot thes pace Bsr} vl 

4 = A, ei elie __ By en i@lglen 

und 

7 A; ei @lales ct B, e 1 Mileg = Ay ei @laler uN B. e t@liler : 

— 

. Die Berechnung fangt beim Intervall 1 an, wo A, = B, fiir x = 0 ist. Man 


wahlt ein passendes w als die gesuchte Eigenkreisfrequenz und nimmt einen 


beliebigen Wert fiir die Vektoren A, = B, an. Bis x =1, gehen sie iiber in 


SULZER 
74541 


Figur 19 
Das (2 — 1)-Massen-System stellt eine Nockenwelle dar. 


A, eles bzw, B, e~'°4/% (siehe Figur 20). Man bildet am Vektor A, e?®4/* ein 
Eechtwinkliges Wheieck mit tg y= w O; ¢c/J,1 G; seine Hypokenies ist dann 
der Vektor 
‘ 


_ @ 0, Cy SPIE 

1 - + 1) Aj Cae 
‘ ( Jruks - 
Auf gleiche Weise lasst sich 
: . OO, cy iwl,/c, 
4 ge 


konstruieren. Subtrahiert man diesen von jenem und multipliziert ihr Resultat 
mit Cy Jp 1/¢1 Jpg, SO kommt der Vektor 


A, ethics ve B, e twly/Co 


mistande. Die Kontinuitatsgleichung liefert ausserdem den Vektor 
| A, ei elles se B, e iwl,/es : 
la er gleich 
; A et Phles ak B, ete ler 

1 
st. Somit lassen sich A, e*@”/* und B,e~'®"/= aus den beiden zuletzt ge- 
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wonnenen Vektoren ermitteln (vgl. Figur 20). Sie sind in ihren skalaren Grossei 
immer gleich und liegen spiegelbildlich um die reelle Achse, mit ihren Pfeil 
spitzen in demselben Abstand von der imaginaéren Achse wie die Vektore: 
Axe Oa nndibee oe hs 

Da die Vektoren A und B immer symmetrisch in bezug auf die reelle Achs 
auftreten, geniigt es, nur einen von ihnen zu zeichnen, zum Beispiel den Vektor 4 
Somit kann die Darstellung vereinfacht werden. Man zeichnet in Figur 2) 
am Vektor A, e?®"/ das rechtwinklige Dreieck aOb und zieht die waagrecht 
Linie bc und die lotrechte ce durch den Punkt a. 

Wird 


SULZER 
74544-A 


Figur 20 


7 ata en iwly|cs —iwlyce : 
Konstruktion der Vektoren Aye #2 und Bye Ohi ays Jenen A, und B, (A; = By). 


Wol. V, 1954 Torsionsschwingungen 315 


gemacht, dann ist Od = A, e'®"/, Aus diesem geht der Vektor B, e ?@h/% ohne 
weiteres hervor. Auf diese Weise geht man von einem Intervall zum niichsten 


Es ist gewohnlich c; = c,, da die Nockenwelle normalerweise aus einem ein- 
-zigen Material gebaut wird. Ist ausserdem J, ; = Jon, 80 fallt der Punkt d mit c 


b 


SULZER 
745 44-8 


1, 


Figur 21 


Der vereinfachte Vektorplan fiir den in Figur 20 dargestellten. 


susammen. Bei einem Ubergang zweier Leitungsquerschnitte ohne dazwischen- 
liegende Masse ist y = 0; dann schrumpft das Dreieck abc zum einzigen Punkt a 
usammen. Der Punkt d wird einfach durch das Verhaltnis de/ae = cy J,4/¢, Jno 
estimmt. 

Figur 22 zeigt den ganzen Vektorzug vom Intervall 1 bis 7. Wenn w richtig 
rraten ist, fallt der letzte Vektor A, e'®!”/ in die horizontale Lage. Ist dies nicht 
er Fall, so muss man es mit einem anderen w versuchen. Mit etwas Ubung 
kommt man auf diese Weise schnell zum Ziel. 


a 


ee ee ee ee 


Die, el tee eek B® bh 


SULZER 
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i al, 
Ane’ &n 


Figur 22 


Der ganze Vektorzug vom Intervall 1 bis fiir das System der Figur 19, 


te ee Be, 
. 
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9. Zusammenfassung 


Die partielle Differentialgleichung fiir ein Torsionsschwingungssystem, unter 
Beriicksichtigung sowohl der Masse wie auch der Dampfung der elastischen 
Glieder, lasst sich streng lésen, wobei die Schwingungen als hin- und hergehende 
Wanderwellen behandelt werden. Die Schwingungen lassen sich wegen der 
Wirkung der Dampfung nicht durch reine Kreisfunktionen ausdriicken. 

Die erzwungenen Schwingungen mit Erregungen an mehr als einer Schwung- 
masse kénnen nur durch das Charakteristikenverfahren gelést werden. Ein 
Mehrmassensystem, wie eine Schiffsdieselanlage,.lasst sich in zwei Teile zer- 
legen. In einem Teil sind erregende Krafte vorhanden, im anderen nicht. Diese 
Trennung gestattet eine wesentliche Vereinfachung der ganzen Schwingungs- 
berechnung. 

Fir Systeme ohne Dampfung wird ein graphisches Verfahren zur Bestim- 
mung der Eigenfrequenz entwickelt, welches einfacher als die bisher bekannten 
rechnerischen Methoden ist. 


Résumé 


Un systéme mécanique composé de masses discrétes reliées entre elles par des 
arbres élastiques dont ni la masse ni l’amortissement ne sont négligeables, est 
susceptible de vibrations de torsion. Dans le cas le plus général, le probléme se 
résoud — d’aprés D’ALEMBERT — par deux ondes élastiques opposées dans leur 
sens de propagation. Du fait de la présence d’un amortissement dans les arbres 
élastiques, les vibrations ne se laissent plus exprimer par de simples fonctions 
goniométriques. 

Les vibrations forcées entretenues en plusieurs endroits du systéme se ré- 
solvent par la méthode des caractéristiques. Tout systéme & masses multiples, tel 
que les installations motrices des bateaux, se laisse scinder en deux parties, l’une 
comprenant toutes les perturbations, l’autre n’en contenant aucune. Cette dis- 
jonction améne une simplification appréciable du calcul des vibrations. 

Un procédé graphique de calcul de la fréquence propre a été établi pour tous 
systémes 4 amortissement nul, mais od la masse des arbres élastiques n’est pas 
négligeable. Ce procédé est plus expéditif que toutes les méthodes analytiques 
connues. 


(Eingegangen: 17. September 1953.) 
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; On the Dynamics of Plastic Circular Plates’) 


By H. GEOFFREY HoPxKINs and WILLIAM PRAGER, Providence, R.1., Wes AL) 


1. Introduction 


The general problem of plastic deformations in structures subjected to 
dynamic loading is very complex. So far only the simplest possible physical 
‘situations have been treated theoretically. These include the propagation of 

plane waves along a rod ([1] to [4])%), the propagation of spherical waves in an 
‘infinite medium [5], and transverse impact on cables and membranes [6], [7], 
beams ([8] to [10]), and frames [11]. Since an exhaustive bibliography of this 
work is beyond the scope of this paper, only a few representative references 
are given. 
_ The only theoretical work on transverse impact on plastic plates of which 
the authors are aware is that of BAKHSHIYAN [12] and Kocuerkov [13]. These 
authors seem to be primarily interested in the piercing of armour plate by a pro- 
jectile. They neglect all bending stresses in the plate, taking account only of 
‘the shearing stresses along the generators of coaxial cylinders around the path 
of the projectile which is supposed to be normal to the plate. In the present 
‘paper, the usual assumption of the theory of thin plates is adopted according to 
which the bending stresses and shearing stresses parallel to the surface of the 
‘plate far exceed the shearing stresses normal to the plate surface. Although the 
‘shear force resulting from the latter shearing stresses is taken into account in 
the equation of dynamic equilibrium, these shearing stresses are neglected in 
the yield condition. 
_ Ina previous paper [14], the authors discussed the load carrying capacities 
of thin circular plates made of a plastic-rigid material that obeys TRESCa’s 
yield condition and the associated flow rule. Various conditions of rotationally 
symmetric loading and edge support were considered. In the present paper, 
the same type of analysis is applied to discuss the dynamic behavior of a 
simply supported circular plate subjected to a uniformly distributed load which 
is brought on suddenly and, after being kept constant during a certain interval 
of time, is removed suddenly. 


1) The # cbaalis presented in this paper were obtained in the course of research sponsored by the 
Ballistic Research Laboratories, Aberdeen Proving Ground, under Contract DA-19-020-ORD-798. 

2) Brown University. 

3) Numbers in square brackets refer to the Bibliography, page 330. 
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2. Equation of Dynamic Equilibrium. Yield Condition and Flow Rule 


Let 7, y, z be cylindrical coordinates, the origin being at the center of the 
undeformed plate’s middle surface, and the z-axis pointing vertically down- 
wards. Furthermore, let R be the plate radius, w the deflection, M the bending | 
moment caused by the radial stresses, and N the bending moment caused by | 
the circumferential stresses!). The variation of the rotationally symmetric 
states of stress and deflection will then be specified by giving w, M, and N as” 
functions of the radial coordinate 7 and the time ¢. 

The equation of dynamic equilibrium of the plate is readily obtained simply 
by including the inertia load in the equation of equilibrium given in the earlier 
paper [14]. If the mass per unit area of the middle surface is denoted by mw, the 
inertia load equals —y w,,, and the equation of dynamic equilibrium takes the 
form (see equation (6) of [14}): 

((¥M),-N=rQ= le [eb W,,) 7 ar, (1) 


0 


where # is the applied load intensity and Q the shear force transmitted across 
unit arc length of the cylinder of radius 7. 

Equation (1) is independent of the mechanical properties of the plate 
material. As in [14] these properties are specified by the yield condition and 
flow rule. The yield condition specifies those states of stress M, N under which 
plastic flow can occur. For each of these ‘yield states of stress’, the flow rule 
makes a statement regarding the rates of curvature 


4=—W,p,,, A= Wert (2) 


that can occur under the considered yield state. 

The following geometrical representation of yield condition and flow rule is 
convenient. Let the yield state M,N be represented by the point with the 
rectangular Cartesian coordinates M, N in ‘stress space’. The locus of the 
points representing all possible yield states is called the ‘yield locus’; it repre- 
sents the yield condition. If x, A are rates of curvature that are possible under 
the yield state M, N, associate the ‘flow vector’ with the components x, A 
along the coordinate axes in stress space with the point M,N on the yield 
locus. The flow rule then establishes a geometrical relation between the points 
of the yield locus and the associated flow vectors. 

The yield condition and flow rule used in this paper are obtained from 
TrEsCa’s yield condition and the associated flow rule. For details of the deri- 
vation the reader is referred to [14]; the results are as follows. 


oie Stee subscripts will be used in the following to indicate differentiation, the notation M,N 


seems preferable to the usual notation M,, Mg. Similarly, x, A will be used instead of Hy, Hq to 
denote the rates of curvature. 
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The yield locus is the hexagon ABCDEF in Figure 1 where M, denotes the 
bending moment necessary to bend the plastic plate into a cylindrical surface. 
If the considered yield state is represented by a point P other than a vertex of 
the yield hexagon, the flow vector is normal to the side containing P and is 
‘directed towards the exterior of the hexagon. If the considered yield state is 
‘Tepresented by a vertex Q of the yield hexagon, the flow vector may have the 
direction of the exterior normal to either of the two sides meeting at Q or to 
any direction between these two limiting directions. 


Figure 1 


Tresca yield hexagon. 


It will appear later that only yield states represented by points on AB will 
occur in the problem discussed in this paper. Accordingly, only the following 
regimes need to be considered (Figure 1): 


Regime A: MeN=Me 020,42 0; (3) 
Regime AB: 0< M< My, N=M,; x=0,420; (4) 
Regime B:7: M=0, N=M,; 02x2-A. (5) 


There is also the possibility that part of the plate is not yielding at the consid- 
ered instant. Since elastic deformations are neglected in the present analysis, 
such a part moves instantaneously as a rigid body, i.e. the rates of curvature x 
and A vanish. According to (2) this means that, at the considered instant, the 
rate of deflection w, is independent of 7. 


3. Circles Separating Different Plastic Regimes. Hinge Circles 


If the entire plate were at all times in the same regime, the determination 
of M, N, and was functions of 7 and ¢t would not constitute a difficult problem. 
For most problems, however, the plate is divided into a central circular region 
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together with a number of surrounding annular regions, each of these regions 
being in a different regime. Moreover, the number of regions and the position 
of their boundaries may vary with time. It therefore becomes necessary to 
establish relations between the mechanical quantities on the two sides of a 
circle separating two regimes. i 

Let @ be such a circle. Dynamic equilibrium requires that the radial bending | 
moment M and the shear force Q are continuous across (© but the circum= 
ferential bending moment N may be discontinuous. The cohesion of the plate | 
material requires that the deflection w and the velocity w; are continuous” 
across © but the acceleration w;, may be discontinuous as will be shown 
below. Moreover, the rates of curvature x and / may be discontinuous across @ 
If A, and hence w,,, is discontinuous across © this circle is called a ‘hinge 
circle’. In the present analysis hinge circles replace the yield hinges of the 
corresponding analysis of plastic beams and frames. 

A hinge circle may be considered as the limiting case of a very narrow 
annulus across which w,, varies rapidly although in a continuous manner. As 
the width of this annulus tends to zero, the ratio x:4 in the annulus becomes 
infinite. Consequently the regime at a hinge circle must be one of the following: 
A, C, D, F, CD, or FA. It follows that w,, must be continuous except, at most, 
at circles where |M| = M,. 

The preceeding discussion shows that in principle a hinge circle need not be 
a regime boundary and conversely a regime boundary need not be a hinge 
circle. However, for the example that is to be discussed in the following, it 
turns out that except at the plate edge all regime boundaries are in fact hinge 
circles and vice versa. On account of this we shall not trouble to distinguish 
between regime boundaries and hinge circles in the following. 

In general the radius @ of a hinge circle © varies with time; the hinge circle 
is called ‘stationary’ whenever 9, = 0. To discuss possible discontinuities of 
the derivatives of w across a hinge circle () represent the variation of o by the 
curve J” with the equation 7 = o(t) with respect to rectangular Cartesian coor- 
dinates 7, t. Since the deflection w must be continuous across J”, the derivative 
of w along J’is also continuous. Thus, 


[w,] i o,|@,| me! (6) 


where the brackets indicate the discontinuities of the enclosed quantities across 
the curve J”. Now w, is continuous across J"; equation (6) therefore shows that 
the slope w, can be discontinuous only across a stationary hinge circle. If the 
hinge circle @is not stationary, w, must be continuous across J”; hence 


[w,:]+ 0{v,,]=0, (7) 
where by hypothesis 9, + 0. 


— 
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: Let us now return to the general case of a stationary or nonstationary hinge 
circle. From the continuity of the velocity w, across J’ it follows that 


[@2] + ere] = 0. (8) 


Thus, the acceleration w,, must be continuous across a stationary hinge circle 
but discontinuous across a nonstationary hinge circle because 2, and hence w,,, 
is by definition discontinuous across such a circle. 

Dynamic equilibrium requires M to be continuous across J’; thus 


[M,] + e:[M,]=0. (9) 


In deriving equations (6) to (9), we have tacitly assumed that ( is the only 
hinge circle that occurs in the considered portion of the plate. If we had to 
deal with two hinge circles Cand (@’, one moving faster than the other and 
eventually catching up with it, we should have to give special consideration to 
the conditions at the point of intersection of the curves J°and I” that represent 
‘the motions of the two hinge circles. 


4. The Mathematical Problem 


a) ey 


We are now in a position to formulate the mathematical problem of the 
motion of a circular plate under uniformly distributed, time-dependent loading 
= p(t). It is assumed that the plate is initially flat [w(r, 0) = 0] and at rest 
{w,(r, 0) = 0], and the origin of time is so chosen that the plate begins to yield 
at t= 0. The deflection w(r, t) and the bending moments M(r, ¢) and N(r, ¢) 
are then to be determined in such a manner that certain conditions are fulfilled 
throughout 0 <7 < R forallt = 0. In enumerating these conditions, we shall 
first assume that (t) is continuous. We then have the following conditions. 

(i) w and w, are continuous in 7 and ¢; the derivatives W,, Wy, Wyy, Wey are 
piecewise continuous in 7 and continuous in ¢, y-discontinuities occurring 
only at hinge circles. Across a stationary hinge circle w, and hence w,, and 
w,, may be discontinuous, but w,, must be continuous. Across a non- 
stationary hinge circle w, must be continuous, but w,,, @,,, and w,, may 
have discontinuities satisfying the equations (7) and (8). 

(ii) M is continuous in and t; the derivatives M, and M, are piecewise con- 
tinuous in 7 and ¢, r-discontinuities occurring only at hinge circles. These 
discontinuities must satisfy equation (9). N is piecewise continuous in 7 
and continuous in ?. 

(iii) w, M, and N satisfy the equation of dynamic equilibrium, equation (1). 
At the center of the plate M = N = M, unless a neighborhood of the 
center moves as a rigid body with w,, = 0 identically in 7. 
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The point M, N lies inside or on the yield hexagon of Figure 1. In the first 
case w,, — 0 throughout the region in which the state of stress is repre- 
sented by interior points of the yield hexagon. In the second case, the rates | 
of curvature and the bending moments must satisfy the flow rule [equations | 
(3) to (5)]. | 
(v) At the simply supported edge, M(R, t) = 0, w(R, t) = 0 for all ¢ 2 0. | 
(vi) Initially, w(7, 0) = 0, w,(7,0) =OforO S7rSR. ; 

When /(t) is allowed to have discontinuities, all quantities considered above, 
except w and w,, must be allowed to have corresponding discontinuities in ¢. : 


= 


(iv 
| 
| 
: 


5. Simply Supported Plate under Uniformly Distributed, Rectangular 
Load Pulse 


Consider a simply supported plate under the action of a uniformly distrib- 
uted load that is brought on suddenly at ¢ = 0, maintained at the constant 
intensity p during 0 < ¢ < t, and suddenly removed at ¢ = Tt. It follows from 
the limit analysis presented in the earlier paper [14], that dynamic plastic 


action will occur only if 
M, 


p> by = oa? (10) 


for p = py the plate yields indefinitely slowly and inertia forces do not arise. 
When (10) is satisfied, the dynamic behavior of the plate is found to have 
different character according as the load is ‘medium’ (f) < p< 24,) or 
‘high’ (pb > 2 pp). 


a) Medium load (py < p < 2 py) 


Two phases must be distinguished in which the motion of the plate is 
governed by different equations. 

Phase I (0 <¢ <1). In the indefinitely slow yielding under the load f, the 
originally plane middle surface of the plate is transformed into a cone (see [14}). 
Guided by this fact consider the motion described by 


w(r,t) = Wie) (L— p), (11) 


where W(t) is a function of ¢ as yet undetermined. This motion satisfies the 
boundary condition w(R, t) = 0 [see (v) above], and the initial conditions (v1) 


will be satisfied if 
WO) = (0) = 0° (12) 
According to (2) the rates of curvature associated with (11) are 


x=0, hawt. (13) 
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When W, > 0 the flow rules (3) to (5) therefore require that 
0=<xM=<M,, N=M,. (14) 


With this value of N and the condition that M = M, for r = 0, the equation 
of dynamic equilibrium (1) furnishes 


M=M,-—~pr++ Wer eae (15) 


At the simply supported edge M(R, ¢) = 0 [see (v)]. This condition and the 
initial conditions (12) will be fulfilled if : 


W(t) = f= fe. £2 (16) 


Note that this gives W, > 0 as was assumed above. Substitution of (16) into 


(15) yields 7, = " " : 

) nite e (Peg). a 
As r varies from 0 to R, the bending moment (17) is found to decrease mono- 
tonically from M, to 0 if f) << p< 2 py. The first relation (14) is therefore 
satisfied, and w as given by (11) and (16), N = M,, and M as given by (17) 
constitute a solution of the problem for 0 < ¢ < rand fy < p < 2 yp. 

| Phase II (rx <t< T). At t=rT the load # is suddenly removed, but the 
plate will not cease moving until the kinetic energy acquired under the action 
of p has been dissipated in plastic work. The instant at which the plate comes 
to rest will be denoted as t= 7. During this second phase of the motion, 
Bb — 0. 
- It so happens that all conditions of the problem can be satisfied by assum- 
ing w to have the form (11) and M to have the form (15) with p= 0. The 
function W(t) for the second phase is determined from the conditions that w 
and w, are continuous at ¢ = t and that M(R, t) = 0 fort < ¢ < T. Thus 


art ciel ae Piss 
Wit) = Fi t(2¢—t) z eee (18) 
M v2 vy 
Mea Re (2- R)- m) 


The expression (19) satisfies the first relation (14) for0 Sr S R, and equation 
18) shows that W, > 0 only when 


Ree kc 20 
oT aa (20) 


At the instant ¢=T the entire plate comes to rest. The final permanent 
leflection of the plate is found by replacing ¢ in (18) by T as given in (20) and 


~ te wry 
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substituting the result into (11). Thus 


w(r, T) = Pip be 7 (1~ 7). (21) 


In particular the central permanent deflection is given by 
w(0, 7) =——, (22) 


where J = pris the total pressure impulse acting on the plate, and [* = (p — fp) t 
may be called the ‘excess pressure impulse’, the term “excess pressure’ referring 
to the excess of the pressure # over the minimum pressure pf, required to 
produce plastic deformation. 


b) High load (p > 2 fo) 
Equation (17) yields 


/] 
M,(0,1)=0, M,,(0,4) =2 8 (2 —2). (23) 
These relations show that the value M(0, t) = M,isamaximum onlyif p < 2 po, 
but becomes a minimum when p > 2 9. In the latter case, the bending moment 
M would assume values in excess of M, in the neighborhood of 7 = 0. Since 
this would violate the yield condition, the expression (17) no longer represents 
a solution when # > 2 f. Since M has the tendency of increasing beyond My 
near the center of the plate, but is prevented from doing this by the yield con- 
dition, we must expect that M = M, over a central circular region 0 <7 < Qo 
when p > po. Thus the regime A applies for the part 0 <7r< gp, and the 
regime AB for the rest of the plate. Since M = N = M, for the regime A, the 
equation of dynamic equilibrium (1) shows that 


m= =, One os (24) 


All points of this central region thus move with an acceleration which is inde- 
pendent of 7 and ¢. 


In the present case the motion of the plate consists of three phases. 


Phase I (0 <t< rt). The preceding discussion suggests that in this phase 
the deflection is of the form 


Pe, 4 (25) 


SONS SE Te 
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where Q, is independent of ¢ but depends on # in a manner that has yet to be 
‘determined. The deflection (25) satisfies the boundary condition w(R, t) = 0 
and also the initial conditions w(r, 0) = w,(r,0) =0 provided that W(t) is 
‘chosen so as to satisfy (12). 

The rates of curvature associated with (25) are 


i x= A= 0). OS7 <4; | 
A as W,(t) (26) 
x= 0, date yy? <7 SR. | 


‘Thus A is discontinuous at 7 = Qy, and the circle 7 = gy is a stationary hinge 
circle. Equations (26) are compatible with the assumption that the regime A 
applies inside, and the regime AB outside, this circle, except for the boundary 


circle y = R which is in the regime B because M(R, t) = 0. We therefore have 

4 M=N=M,, 0=r<o, 0<XM<™M,, | 

| : (27) 
W=M,, o<*<R, M=0, N=M,, r=R. | 


From (24), the first equation (25), and the initial conditions (12), it follows that 
A 2 
j wiy=2F. (28) 


The bending moment M must be continuous across the hinge circle 7 = @p. 
Thus M (0p, t) = My; moreover, M(R, t) = 0. With these boundary conditions 
and with w(r, t) as specified by (28) and the second equation (25), the equation 
‘of dynamic equilibrium (1) applied to 99 [7 S R yields 


oS ee (29) 
and 
M ras | j ea he (vy + Qo) os pe he (30) 


According to (30), M decreases monotonically from M, to 0 as 7 increases 
from g, to R. The continued inequality for M in (27) is therefore fulfilled. 

Equation (29) determines @)/R as a function of p/(2 po). The value of @9/R 
corresponding to a given value of p/(2 f)) can be read off Figure 2. As # in- 
creases indefinitely, 09 is seen to tend towards R. 

Phase II (rx < t < T,). During this phase the radius @ = g(t) of the hinge 
circle is found to decrease from @(t) = @9 to e(Z,) = 0. Indeed if we assume 
that the regime A continues to apply near the center of the plate the reasoning 
which led to equation (28) shows that the diminishing central portion of the 
plate must preserve whatever velocity it had at ¢ = t. If the radius of the hinge 
circle remained constant under these circumstances, the kinetic energy of the 
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Figure 2 


Variation of p/(2 po) or T,/t (ordinate) with @9/R or o(t)/R (abscissa), respectively. 


plate would remain constant in contradiction of the fact that energy is dissi- 
pated in plastic deformation. In view of (28), these remarks suggest that the 
velocity during the second phase is of the form 


| Puc R—y Cu 


uw  R— eit)’ a(t) 


where g(t) remains to be determined. The reasoning which led to the equations 
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ST aE eee ae 


(29) and (30) now yields the following relations: 


mae 

: Q(B +2 Rp 3%) = SePoR (32) 
. pt 

and 

. 

Mu Pee) Pane el ree oR 3 Oi 

4 “ao ie anes v (R — 0)® (R + 30) pO) 7 Se 33) 
, ) 

_ __ Equation (32) must be integrated under the initial condition that o(t) = gp. 


With the use of (29), we thus obtain the following cubic equation for o = o(t): 


4 He) 3 @ 2 Lea 2Pot 
4 (ze) =(z) CT ty hae a (34) 
Since o(7,) = 0, by the definition of 7,, equation (34) yields 
a 
‘ pt 
T. —— 
eG ge Fo) 


yer ee 


a 
_ and the right-hand side of (34) can be written as (¢/T,) — 1. It is then seen that 
- equation (34) establishes exactly the same relation between 7,/t and 0/R as 
ee equation (29) between #/(2 #9) and @,/R. Thus a simple relabeling of the 
_ coordinate axes in Figure 2 will enable us to find, from this figure, the value 
4 of o/R that corresponds to a given value of T;/t. 

_ For a fixed value of t, and hence of @, the bending moment (33) decreases 
monotonically from M(o, t) = M, to M(R, t) = 0. The continued inequality for 
_M in (27) is therefore satisfied. 

. Finally the deflection w(r, ¢) remains to be determined. We shall here be 
concerned only with the deflection w(r, T,) at the end of the second phase, and 
4 we shall determine this by making use of the jump condition (7). 

. Equations (25) and (28) show that, at the end of the first phase, the middle 
_ surface of the plate has the shape of a truncated cone, the deflection of the points 
inside the hinge circle y = @ being # t?/(2 w). At the end of the first phase, a 
discontinuity in slope has therefore developed at 7 = Qp: 


; Bie 
W,(Oo a O57) = (Oo 0st) Die 2 Gai (36) 


During the second phase the hinge circle is contracting. Since the velocity 
-w, must be continuous across this nonstationary hinge circle, equation (6) 
- shows that new discontinuities in slope cannot develop during the second phase. 
According to the first equation (31), the contracting hinge circle moves into 
a region in which w,, = 0. Consider a fixed circle ( of the radius ry. As the 
‘contracting hinge circle moves across ( it produces there a curvature w,, 
_of the plate meridian which can be found from equations (7) and (31). Indeed, 
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according to (31), we have at H 


z ) 
[ers] a ae ay (37) 

and hence, by (7) and (32), 
[w,,] ca f= £ ee ss — : (38) 


Since GH moves into a region in which the plate has performed a vertical 
translation only, w,, vanishes inside A The jump (38) thus equals the value 
of w,, just outside WH The second equation (31) shows that w,,; = 0 outside Hy 
so that no further change in the curvature w,, occurs at C once. has moved 
across @ Thus, 

(emt) hee Shy 


pL ee (39) 


DyplT, t) = 2 ‘ Po : R3 


At the end of the second phase, equation (39) applies right up to the center 
of the plate. Since 


2 
w,(0,T)=0, w(0, T) = 5 Es 1), (40) 


by (28) and the first equation (31), integration of (39) with the initial conditions 
(40) yields 


ie T eee ES (ae ntaes ie —1], 0<r<q. (41) 


Outside the circle of radius 79, the middle surface of the plate assumes the shape 
of a right circular cone. This is determined by the deflection w(o9, T,) found 
from (41) and the condition that w(R, T,) = 0. We leave it to the reader to 
check that the discontinuity in the slope w, of the middle surface at 7 = ay 
when ¢ = T, has the value (36). 

Phase IIT (T, <t < Tj). At the end of the second phase the entire plate 
except the center and the edge is in the regime AB. The subsequent motion 
of the plate is therefore described by equation (11) in which W(t) must be 
suitably determined. The analysis of this third phase of the motion proceeds 
along similar lines as the analysis of the second phase in the case of medium 
load. We leave it to the reader to work out the details and state only the 
following facts. 

The plate comes to rest at the time 


T,=27,. (42) 


=~ 
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The deflection at this instant inside the circle of radius Qo 1S given by 


. PET <p y ae 73 
Oe T= ag. 3p ee ar) 1) 6 OST SO, (43) 


and outside this circle is linear in r (09, <r < R), vanishing of course at the 


Tete ee eae 


aight oy 


; 


+ helt ett 


_ plate edge. In particular it follows that the final central deflection is 

4 
lage i 
4 0, Z;) = —_,—-"—__, 

, oleae 

4 


where I = # tis the total impulse acting on the plate. 
j It should be noted that the meridian of the permanently deformed plate 
is curved for 0 < 7 < g, and straight for gy < 7 < R with nonvanishing slope 
at y = 0 and a discontinuity of slope at 7 = Qo. 

In Figure 3, the variation of w(r, T,)/w(0, Tz) with 7/R is shown for some 
-Tepresentative values of f/(2 9); as #/(2 Ho) increases the central curved part 
of the permanently deformed plate is seen to increase in size. 
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Figure 3 
Variation of final deflection across plate, expressed in terms of central deflection, with pressure 
| intensity. (The circles indicate the positions of the hinge circles in phase I.) 


Abstract 


The paper is concerned with the dynamics of a thin circular plate made of a 
plastic-rigid material that obeys the yield condition of TREsca and the associated 
flow rule. The basic equations are established in the case of rotational symmetry, 
and their complete solution is given for a simply supported plate subjected to a 
uniformly distributed load which is brought on suddenly and, after a certain time 
interval, removed suddenly. 
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Résumé 


Cet article étudie la dynamique d’une plaque circulaire infiniment mince 
fabriquée avec un matériau plastique qui obéit a la condition de plasticité de 
TREscA et a la loi d’écoulement associée. Les équations fondamentales sont 
établies dans le cas ot existe une symétrie de révolution. Leur solution complete 
est donnée pour une plaque reposant sur des appuis simples, soumise a une charge 
uniformément répartie qui est appliquée rapidement et demeure constante 
pendant un certain temps, apres lequel elle est rapidement retirée. 


(Received: October 2, 1953.) 


Stress Distributions in Semi-Infinite Solids of Revolution 


By Antony J.A. MorGan!), Pasadena, California 


Introduction 


All of the problems presented in this paper are solved by reducing the govern- 
ing partial differential equations to ordinary differential equations for which 
closed solutions are found. This reduction is accomplished by the application of 
a recently developed theory [1]*) for obtaining similarity solutions’) of systems 
of partial differential equations. The term invariant solutions is suggested in 
place of similarity solutions because this terminology reveals their underlying 
mathematical nature [1], [2], [3]. 

The solution to the problem of obtaining the stress distribution in a solid 
bounded by two concentric cones with a common apex, under external loadings 


1) Aercon, Inc. 
2) Numbers in brackets refer to the References, page 340. 
8) This terminology is used in aerodynamics (laminar boundary layer, conical flow, etc.). 
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which vary as an arbitrary real number power of the distance from the apex, 

has not been previously reported in the literature. Solutions of this problem 

are presented in this paper. 

If « > B are the semi-vertex angles of the two bounding cones, three cases 

are distinguished: 

(i) « > a/2 > B: this can be considered as a circular plate of infinite radius ; 

(ii) 2/2 > « > B: corresponds to a semi-infinite conical shell, and 

(iii) 2/2 > « > 0, one bounding cone: is equivalent to a semi-infinite solid 

cone. 

In all of these cases the thickness increases linearly from the apex. These 

E distinctions, (i), Gi), and (iii) (Figures 1, 2, 3), are of importance since they 
“dictate different choices for the form of the stress-function solution. 
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Nomenclature 


The following notation is used in the paper: 


a numerical parameter in a group of transformations; 

BC,,---,C, constants of integration; 

D differential operator = d/dy if y is the independent variable; 

d greatest value of 7 or z, at which one of the proportional limits 

of the material is attained; 

e least value of y or z, at which one of the proportional limits of the 

, material is attained ; 

LF (6) dependent variable of equation (15), defined by equation (12) ; 

-F(o@,6;7; 2) hypergeometric function of real argument z; 

—G(n) new dependent variable for equation (15), defined by equation (22) ; 

R tan «; 

‘R’ tan B; 

mM, Nn, numerical constants associated with the group in which a is the 
parameter ; 

t_ normal loading (force/unit area) on upper face of semi-inf.nite 
solid of revolution; 

1c) Legendre function of the first kind of order w and argument ¢; 

q a real number exponent denoting the type of variation of the 
external normal and tangential loads; 

Q,() Legendre function of the second kind of order w and argument ¢; 

er tangential (shear) loading (force/unit area) on upper surface of 

; semi-infinite solid of revolution ; 

v(C) new dependent variable for equation (15), defined by equation (I. 2) ; 

7,0,2 cylindrical coordinates ; 

4 angle of upper surface of semi-infinite solid of revolution with 


respect to r- or z-axes (Figures 1, 2 Oks 


O,, O69, Oz, Trz 
P(r, z) 


(a0) 
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angle of lower surface of semi-infinite solid of revolution with 
respect to 7- or z-axes respectively (Figures 1, 2) ; 

q+ 3; 

the gamma function; 

& (1 + £%)-4/, for circular plates; 

(1 + 7?)~¥?, for conical shells and solid cones; 

(r/z) = similarity variable used in semi-infinite solid of revolution 
(conical shell) problem ; 

q+1; 

Poisson’s constant; 

(z/7) = second similarity variable used in semi-infinite solid of 
revolution (circular plate) problem; 

stresses, cylindrical coordinates (axial symmetry) ; 

stress function in cylindrical coordinates (axial symmetry) ; 
angle of ray within semi-infinite solid of revolution with respect 
to r- or z-axes (Figures 1, 2, 3). 


Primed quantities denote external loadings on lower surface of semi-infinite 
solid of revolution. 

Dashes denote differentiation with respect to independent variable indi- 
cated by the context. 


Circular Plates of Infinite Radius and Linearly Varying Thickness 


This corresponds to case (i) (Introduction), an axially symmetric problem, 
where the external loadings are as shown in Figure 1. 


Tetse 


Figure 1 


The Circular Plate of Infinite Radius and Linearly Varying Thickness. 
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The stress distribution within the plate is determined by a solution of the 
partial differential equation [4], [5] 


0 cali 0 02 


a2 
VP = (sat 3 or T on a +5 “aap hg PPD (1) 


” 


where @ is the stress function in cylindrical coordinates and the stresses are 
given by: 


at the lower face (w = f) of the plate. Note that the angle w is measured with 
tespect to the 7-axis. 

The form of the solution is deduced as follows: Equation (1) is unchanged 
(constant conformally invariant) under the transformations’) 


; 5 ey 
; y7=a"r, Z=a"z, D=a°O, (10) 


- 0 02D 

4 obs [> ear , (2) 
4 f) 1 0 

: ig Oe iaheeg orale (3) 
+ 

~ 0 02D 

; Og = [(2 — v) 72D — YE : (4) 
and 

: 0 02D 

: Tre= Gy [(L= 9) 199 — Ser] (0°) 
Z 

: The desired solution of equation (1) must satisfy the boundary conditions: 
4 P—Th=—o~ toh: (6) 
: Pk+T=—o,k—-%, (7) 
3 at the upper face (w = —a), and 

; PY — Tk’ = —0, 4 beh’, (8) 
PUREST ah Eg, (9) 
; 


‘where a (+ 0) is a numerical parameter and m,n, p are real constants, if 
(n/m) = 1 and (p/m) is arbitrary. Two absolute invariants of the transforma- 


tions (10) are 


@ 
i son and pm * (11) 
| Hence, the substitution of 
P= rmFO, E=2 (12) 


1) Mathematically, this is a continuous one-parameter group of transformations. See [1], [2], 
‘and [3] for terminology used in group theory. 


y 
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into equation (1) will, by a known theorem [1], reduce it to an ordinary differ- 
ential equation with F(é) and & as dependent and independent variables 


respectively. 
Equation (12) gives the desired form of the solution, if the external loadings 
are of the form: 


Pa hr, Nes ES Pe ee (Ey Fig = (13) 
when 
Pee 
a ae (14) 


This is seen on substituting equations (10) in equations (2) to (5) and then 
substituting the results, together with equations (13), in equations (6) to (9). 

The substitution of equations (12) into equation (1) gives the ordinary 
differential equation 


(Pye) 22 gs tcc? ft ee) 

Pl24g-F 1) eth =e) — [2+ @+2)— 12ers eee 
+ (g +1)? (q+ 3)*F() =0 

whose general solution (Appendix I) is: 


F(é) = (1 = £*)- 9) [C, P(t) + €C,0,(0) + C, B() + COC), (ie 


where 
Gomie (Lie ot) Os per iy, 
eg ara 
v= g = 3, 


P,(¢) = Legendre Function of the first kind of order w with argument ¢, and 
Q,,(¢) = Legendre Function of the second kind of order w with argument €. 


By direct calculation, on substituting equations (12) and (14) into equations 


(2) to (5), the stresses, in terms of g, €, and F(€), are given by the following 
expressions: 


1" oy = (q + 2)» (¢ + 2) — (G+) F'(é) 
+ [2 (¢+1) — 9(2q+3)]EF"(E) +[y (1+ &) — 2] FS, 
ra, = (9 + 2)[y (g + 2) — 1] P(e) — 9 (29 + 3) EFC) . 
VR Ree Sieh Pall (ai 


4 
Vol. V, 1954 Stress Distributions in Semi-Infinite Solids of Revolution ; 335 


= (2—») (g + 22°F) — (2-9) (29 + 3) €F"(é) 


Eee ON, 
2 
A 


+([(2 — ») (1+ €) — 1] F"() MW) 
“and 
a =(1— ») (¢ +1) (¢ + 3)? P(E) — (1— ») 3g? + 11¢+ 9) EF 6) 
4 +[(1—») (3944) &—»(¢+ 1] FS) (20) 
E +[1— (1—») €? +1] EF") 
| Equations (13) and (17) to (20) can be used to express the boundary con- 


itions (6) to (9) in terms of Py, Ty, Po, Ty, k, k’, F(€) and the derivatives of 
Fé). The stress distribution within the circular plate (Figure 1), under external 
loadings given by equation (13), is then completely determined by the solution 
(16) of equation (15). 
__ Several cases must be considered: 
4 Case A, q > 0. Here the stresses, equations (17) to (20), become infinite as 
-¥ co for any fixed € within the body. The proportional limits of the material 
(in shear or tension) are then attained at some values of r > 0. Let r = e be 
the least such 7. The problem will retain physical significance if the plate is 
cut crosswise at some 7 < e and the stresses given by equations (17) to (20) 
‘are applied at the cut as external forces. 
_ Case B, q< 0. The stresses are infinite at the origin (ry = 0). Let ry = d be the 
‘greatest value of 7 at which one of the proportional limits of the material is 
‘attained. Physical significance is retained on cutting a hole, of radius > d, in 
the center of the plate and applying an appropriate system of stresses [given 
by equations (17) to (20)] at the cut as external forces. 
The case g = 0 (uniform loading) does not yield infinite stresses at any 
‘point in the plate. The restrictions outlined under cases A and B will also apply 
in all subsequent considerations. 


; 
os 


Semi-Infinite Conical Shells of Linearly Varying Thickness 


This corresponds to case (ii) (Introduction). The situation is illustrated by 
Figure 2. The solution (16) could also be used for this case if the external 
loadings were expressed as a power of 7. For ease of application it is convenient 
to seek a solution which allows the external loadings to be written as proportio- 
nal to a power of z and the cone vertex angles to be measured with respect 
to the z-axis. This is possible since equations (12) and (14) indicate that the 
stress function ®(r, z) may be written as: 


1 
® = 2 F(T), q=—E=p (21) 
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Figure 2 
The Semi-Infinite Conical Shell of Linearly Varying Thickness. 


where F(€) is the solution (16) and y is any real number. Thus the required 
solution is 


G(n) = 1" F (;) ; (22) 
where y = ¢ + 3. Hence, by equation (16), 
G(y) = -? [C, P,(2) + Cz Q,(0) + Cy BS) + C1 0,(2)], (23) 


where, now 4 
Paid) ie phe 2s 


and the other quantities remain as previously defined. 
Since ® = z” G(y) the stresses, equations (2) to (5), can be written in terms 
of G(y), and its derivatives as follows: 


a4o, = ¥(g+1) (g+ 2) (¢ + 3) Gy) — » (¢ + 2)[3 (¢ +1) yn? -1] 9-1 G'(n) 
— [3 9 (9 +1) 2+ gv — (g+1)]G"y) —[v +7) -1]n GE"), 
240, =» (¢ +1) (¢ + 2) (g + 3) Gy) | 
— (¢+ 2) [3 (g +1)? + (1 — »)] 971 G'(m) (25) 
+ [3-9 (941) fg PLY" y) — rp Teme, | 
24a, = (1—») (g +1) + 2) (¢ + 3) Gin) 
—[3 (1 — ») (@ +1) @ + 2) 4? — (2 — »).(g + 2)] 1 G'(n) 
+ [3 (1 — ») @ +1) 4? + (2— ») g] G"(n) 
=| (1 =») 92 (29) GMs 


(24) 


(26) 
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tre = —[ (9 +1) (Gg + 2) 4? + (1 — »)] 9-2 G'(m) 

i 2 " (27) 
+[2»(¢ +1) 9? + (1 — »)]y-2G"(y) + [1 — » (1 + 92)] 6"(y) , 


here the dashes denote differentiation with respect to 4, and q is any real 
number. 
The boundary conditions are: 


P= Th a, 4 te, Bi, (28) 
Pk+T=-o,k+1t,, (29) 
at the outer surface (w = «) and 

Pe Te meget ok, (30) 
P'k'—T'=—0,k'4+1T,, (31) 
at the inner surface (w = B) of the shell. 

Comparison of equations (24) to (27) and (28) to (31) shows that the bound- 


ary conditions can be written entirely in terms of 7, G(n) and the derivatives 
of G(n) if the external loadings are of the form: 


eet et eee Te gt Pt al Poets JT” ear ee (32) 


‘ton the problem can be solved entirely in terms of the solution (23) of equa- 
tion (1). 


The Semi-Infinite Cone 


The problem is that of obtaining the stress distribution within a semi- 
infinite cone (Figure 3). This can be accomplished by placing certain ‘over-all’ 


Figure 3 
The Semi-Infinite Cone. 


; 
: 
j 
-ZAMP V/22 
a 
j 


338 Antony J. A. MorGAN ZAMP 


restrictions on the solution (23). Note that only two boundary conditions, 
equations (28) and (29), are at our disposal, and hence the solution can only 
contain two arbitrary constants. To eliminate the two additional constants 
occurring in equation (23) the physically reasonable ‘over-all’ condition that: 
‘the stresses must be finite on the z-axis’ is imposed. This requires an exami- 
nation of the singularities of the solution (23) as € >1, 4 >0. 

A suitable representation of the Legendre functions, in terms of the hyper- 
geometric function [6], F(0, o mie z), is 


P(t) =F(u +1, ws; — ZS) -(-1<0<3), (33) 
and 

12 Pi +1 Pee 4 a 
Ose) = Seti ea ae 1F(Sut+ sy qetlwt+ sie) (C >1)§ 


(34) 


similar expressions holding for P,(¢) and Q,(¢). The functions Q,, ,(¢) have a 
logarithmic singularity at 7 = 0(¢€= 1) and so will their derivatives with 
respect to [6]. The P,,,(¢) and their derivatives, on the other hand, are 
regular at ¢ = 1. This eliminates the functions Q,, ,(¢) from the solution (23). 

It remains to be verified that the functions P, ,(¢) do not yield infinite 
stresses at the z-axis. Such a situation is possible if 7~1G’(y), contained in 
equations (24) to (27), is singular at 7 = 0, but the expression 


; ae ? 
Miah | ices Sabor a ae 9) 


is not singular as ¢ > 1. Hence the stresses o,, og, o, are regular on the z-axis. 
It can be shown, on expressing the derivatives of G(y) in terms of derivatives 
with respect to ¢, that the terms 7~? G’(y) and 7-1! G"(n) do not make t,,, equa- 
tion (27), singular as 7 > 0. In fact the singularities contributed by these two 
terms cancel out identically. 

With the above verification the desired solution is immediately seen to be: 


Gin) =o" [or F (uw 4, aes é) vs C.F (y +1,-y31;5>5 é)| 


2 2: ae 
(35) 
where, now 
02 S00. ie Gat 


I(t) = Iie +n) P'(o4+n) gn 


BOs Cont 2) Se I(@) I'(c) I(t +n) ““n!\ 


(|z| <1), 


n=0 


and the other quantities remain as previously defined in equations (16), (23) 
and (33). 
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A 
é The problem of determining the stress distribution within a semi-infinite 
‘S solid cone, under external loadings of the type given by equation (32),can then 
a be solved in terms of the solution (35) of equation (1). 

3 . 
> : 
be Conclusions and Suggestions for Further Research 
4 


The solutions of the problems posed can be extended in several ways. By 

“superposition it is possible to obtain solutions for cases where the external 
loadings are expressed as power series in the distance from the apex (Figures 1 
to 3). The stress distributions for solids which are the finite counterparts of 
those shown in Figures 1 to 3 might also be determined, for certain types of 
boundary conditions at the crosswise surfaces, by superposition of the solutions 
given in this paper. Note that the superposition is made by using the stress 
functions and not the solutions of the ordinary (similarity) differential equation 
corresponding to the governing partial differential equation. 
Invariant (similarity) solutions of the non-homogeneous biharmonic equa- 
tion can also be obtained by use of the methods illustrated in this paper [7]. 
As suggested by TRANTER’S [8] success in the case of the semi-infinite wedge, 
it is possible to generalize the problems solved herein to include boundary 
conditions composed of piecewise continuous external loadings. This is done 
by application of the Mellin transform to equation (1), when expressed in 
spherical coordinates. Of course, the loadings are arbitrary to the extent that 
their Mellin transform must exist. The results of such an investigation will be 
' presented at a later date. 


if 
L 


€ 


Sse Pe 


APPENDIX I 
The Solution of Equation (15) 


Equation (15) can be split into the operator factors 


os I ae ee be Seeitntet Ss. eer ae RN: bh ey F 


[(g2 +1) D?+2(¢+1ED+ (GQ +1)" | ca 
x [G2 +1) D8 (29 +5)ED+ (G+ 3)]FO =O, | 
where j 
D= (ie): 


; Two solutions of equation (15) are obtained by setting the second of the 
3 operator factors in equation (I.1) equal to zero and applying the following 
transformations of dependent and independent variables: 


Fé) = (2 +82 (0), ¢=F +E)". (I. 2) 
The new dependent variable v(¢) satisfies the Legendre differential equation [6]: 
(C2 —1) (0) +2000) —y vt) 0) =0 Yagt3), (13) 


| 
| 
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eres | 

which has as two independent solutions the Legendre functions of the first kind, 
P,(é), and of the second kind, Q,(¢)’)?). . 
On applying the transformations (I.2) to equation (15), the following 


operator representation of the transformed equation is obtained: 


[(¢2-1) D?+20D—p(w+1)][(?-1) D?+20D—y(y +1] (t) =0, 
(1.4) 


where 


D= se, p=¢qst1, and y=qt3. 


It is a simple matter to verify that the two operator factors in equation 
(I. 4) are commutative. Now, since each of these operator factors is a Legendre 
differential equation it is immediately seen that the general solution of equation 
(15) is given by equation (16). 


Abstract 


The problem of obtaining the stress distribution within semi-infinite solids 
bounded by two concentric cones with a common apex is formulated in terms of 
a fourth order ordinary differential equation. This is accomplished by using a 
recently developed theory for obtaining similarity solutions of partial differential 
equations. Such a formulation is valid when the external loadings are proportional 
to an arbitrary real number power of the distance from the apex. Closed solutions 
are given for these equations. 

For certain ranges of the cone vertex angles these solids, all of linearly varying 
thickness, can be descriptively classified as: (i) circular plates of infinite radius, 
(ii) semi-infinite conical shells, and (iii) semi-infinite solid cones. 
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Résumé 


Le probleme de déterminer la distribution de la tension dans les solides semi- 
infinis délimités par deux cénes concentriques ayant un sommet commun est 
_exprimé sous forme d’équation différentielle ordinaire du quatri¢éme degré. Ceci 
est fait par l’application d’une théorie récemment énoncée pour l’obtention de 
solutions similaires des équations différentielles partielles. Une telle expression 
est valide quand les forces appliquées extérieurement sont proportionnelles a 
une puissance d’indice réel, choisie arbitrairement, de la distance au sommet. Des 
solutions exactes de ces équations sont données dans le texte. 

Pour certaines limites des angles du sommet d’un céne, ces solides, tous 
d’épaisseur a variation linéaire, peuvent étre classés descriptivement comme suit: 
io plaques circulaires de rayon infini; 2° enveloppes céniques semi-infinies, et 


3° cénes solides semi-infinis. 


; 


a 


ban aad 
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Die Berechnung einiger komplexer Werte auf einer Multiplizier- 
maschine mit nur einem Multiplizierwerk 


Von WALTER BosscH, Ziirich 


I. 


Auf einer Multipliziermaschine mit 3 7 +1 Stellen im Multiplizierwerk und 
2-(3n-+1) Stellen im Resultatwerk konnen Real- und Imaginarteil einer kom- 
plexen Zahl mit ” Stellen durch einen Abstand von 7 + 1 Stellen getrennt werden. 
Addition und Subtraktion sind trivial’). 


1) Regelungstechnik 1953, H. 12, S. 269-292. Zitiert als [1]. 
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Il. Multiplikation komplexer Zahlen 
a) FF, (Ff, = 3, = komplexe Zahl) 


HF = (a, + by 1) (a, + by 1) = 4, a, + (a, by + ag by) 7 — Dy bg. (1) 
Vergleiche (a, + b,) (a, + by) = 4, Gy + (4; by + aq by) + 0, de. (2) 


Durch Koeffizientenvergleich wird ersichtlich, dass man von der Gleichung (2) 
das zweite Glied mit i und das dritte Glied mit —1 multiplizieren muss, um die 
Gleichung (1) zu erhalten. Unter Beriicksichtigung des unter I erwahnten Stellen- 
abstandes der Zahlen ergibt sich auf einfache Weise die Multiplikation komplexer 
Zahlen. 

b) F F® (F und F* = konjugiert komplexe Werte) 


Fi Be = (ayo, 7) (a, — by 2); (3) 
FB, Ff = at — a,b,i+ a,b,14 03. (4) 


Aus [1] geht hervor, dass die Zahlen, die diesem Produkt entsprechen, im Produkt- 
zahlwerk der Multipliziermaschine in der Reihenfolge der Gleichung (4) erscheinen. 
Um unmittelbar die Summe (a? + 6b?) ablesen zu kénnen, wird umgeformt: 


EER = — (a, + 5,2) (b) + a, 1), (5) 
F, Ff = —a, bi + (a2 + bY) + a, d,7. (6) 


Aus [1] geht wiederum hervor, dass die Zahlen auf dem Produktzahlwerk in dieser 
Reihenfolge erscheinen. 
Damit kann aber ganz allgemein auf der Multipliziermaschine die Summe der 
Kathetenquadrate eines rechtwinkligen Dreieckes direkt abgelesen werden. 
Andere Funktionen von ahnlichem Aufbau sind analog berechenbar. 


III. Division komplexer Zahlen 


Fy a+ O38 al + bE 3 a 


fe ao (a, + by 4) (a, + by i) . 
3 e 7 3 4 


Der Nenner wird wie F, F berechnet. Damit ist die Division komplexer Zahlen 
auf eine Multiplikation zuriickgefiihrt. 


» A 
Fy 
Fs helt she, ck F? 
date caine Rarraare (8) 
Fj Varna ay + OF 


Damit ist die Division konjugiert komplexer Zahlen in eine Multiplikation 
(Quadrat) komplexer Zahlen iibergefiihrt worden. 


Ns i i hi 
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‘df ; 
3 IV. Der Kettenbruch mit komplexen Zahlen 
ele! 
setae i (°) 
: F, F, 
jj me [(a + 6, 1) (a, + by ] [(d, + 52) (4, + 42) 4] PY, (10) 
f (4, + a)? + (6, + by) 


4 Der Nenner wird analog F, F% berechnet, damit ist der Kettenbruch auf eine 
-Multiplikation komplexer Zahlen F, = (1/h) F, —, E, vereinfacht. 


V. Praktische Hinweise 


7 


A Die praktische Ausfiihrung einer Multiplikation komplexer Werte ist in [1] 
_beschrieben. Die hier neu angefiihrte Multiplikation komplex konjugierter Werte 
ist nach Studium obiger Literaturstelle ohne weiteres auszufiihren. 


Summary 


It is explained how a normal mechanical calculator can be used in a very time- 
saving manner for all four basic forms of calculation with complex numbers, and 
also a short method of calculating the squares of catheds in a right angled triangle. 


_(Eingegangen: 22. Marz 1954.) 


Prifgerat fiir Photomultiplier 


Von EpGar LUscHeEr, Lausanne?) 


Fiir spektrometrische Anwendungen der Photomultiplier miissen diese zunachst 
gepriift werden, bevor man sie zu Messzwecken brauchen kann. Im folgenden ist 
eine einfache Apparatur beschrieben, die gestattet, den Dunkelstrom und die Emp- 
findlichkeit von Photomultipliern zu messen. Das Gerat wurde vor allem zur ra- 
-schen Priifung von grossen Serien 
von Vervielfacherréhren entwik- 

kelt, zum Zwecke, diese R6hren zu 
klassieren. Die optische Messan- 
ordnung ist in Figur 1 dargestellt. 
Konstantes Licht von der 
Lampe L tritt durch die Matt- 
-scheibe oder den Spalt auf den um 
die Achse Z drehbaren Lichtdif- 
fuser, welcher das Licht auf den 
Vervielfacher streut. Diese An- 
ordnung erlaubt, die totale Emp- 
findlichkeit der Photomultiplier Figur 1 
in Funktion des Drehwinkels «zu —_, Wolframlampe; S Mattscheibe oder beweglicher Spalt; 
bestimmen. Der Lichtdiffuser 4 W Trennwand; M Lichtdiffuser; P Photomultiplier; 
ist mit einer Schicht Magnesium- Z Drehachse der Diffuserscheibe. 


oxyd versehen, um eine méglichst 
gute Lichtzerstreuung zu erhalten. Die Wande des Priifgehauses sind schwarz. 


1) Applied Research Laboratories, Lausanne. 


; 
344 Kurze Mitteilungen — Brief Reports - Communications bréves ZAMP 

Besondere Beachtung muss der Speisung der Lampe L geschenkt werden. Die 
Spannung wird mittels eines Sérensen-Regulators stabilisiert. Zusatzlich sind noch 
zwei Stabilisatoren OA3 geschaltet. Die Intensitat der Lampe kann mit einem 
Widerstand variiert werden. Wird die selektive Empfindlichkeit gemessen, kann 
eine Spektrallampe benutzt werden und an Stelle von S eine entsprechende 
Filteranordnung. Ein Elektrometerverstarker mit dynamischem Kondensator mit 
einem Potentiometerschreiber am Ausgang dient als Messeinheit. Das prinzipielle 
Schema der Messanordnung ist in Figur 2 wiedergegeben. 


TEV 


O © 
kegistrierendes 
papal aly 


O) 
—-0 
Figur 2 


Schema der Messanordnung. 


Mit den Widerstanden R,, R,, R,; kann die Empfindlichkeit des Verstarkers 
eingestellt werden. Zur Konvertierung der Gleichspannung in Wechselspannung 
dient der elektromagnetisch angetriebene Plattenkondensator DC. 

Die Doppeltriode 5691 mit dem Transformator T dient zur Phasengleich- 
richtung. Die Réhre 6SK7 dient zur Riickkopplung. Mit dem Potentiometer R, 
kann der Nullpunkt des registrierenden Potentiometers verschoben und somit 
dem Messbereich angepasst werden. 

Mit einer Eichspannung am Verstarkereingang kann das Gerat geeicht wer- 
den; damit kann nachher die absolute Grésse des Dunkelstromes bei maximaler 
Dynodenspannung gemessen werden. 

Der Bestimmung der totalen Empfindlichkeit muss eine Eichung fiir den 
Drehwinkel « vorausgehen. Der Ausschlag am Potentiometerschreiber wird in 
Funktion des Drehwinkels bestimmt und die Réhren danach klassiert. 


Summary 


Equipment for testing photo-multipliers with diffused light and with a special 
electrometer amplifier with dynamic condenser is described. 


(Eingegangen: 20. Februar 1954.) 
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Asymptotische Entwicklungen fiir Fresnelsche Integrale und verwandte 
Funktionen und ihre AnwendungsmOglichkeiten bei der Berechnung 
spezieller Raketenbahnen 


Von JOSEF ZBORNIK, Chur?) 


, Wie H. HOsLI in einer bemerkenswerten Arbeit zeigte?), lasst sich unter den 
‘itiblichen Annahmen der Windeinfluss auf Raketenbahnen verhdltnismassig ein- 
fach mathematisch erfassen. Die numerische Auswertung allerdings bereitet inso- 
_ fern einige Schwierigkeiten, als die Losungen auch Integrale der Form 


x 
ae 2d (§=V—1,0Sk S1,10 S x S 100) (A) 
0 


a 
me weisen, deren mégliche Bestimmung mittels der tiblichen Potenzreihenent- 
4 wicklungen praktisch versagt, weil man in Anbetracht des gegebenen Bereiches 
von *# eine betrachtliche Anzahl von Gliedern berechnen miisste, um die ge- 
_ wiinschte Genauigkeit?) von + 2- 10-4 zu erreichen. 

Diese Schwierigkeit kann durch Entwicklung des Integrals (A) in eine asym- 
ptotische Reihe vermieden werden, wie die nachstehenden Ausfiihrungen zeigen. 


Das Doppelintegral 


J uty he 

} =V2 f[ [ .-0*+*-D aya (k >0, reell 

s eel Ff fe ydt (k=O, reel), 

4 i) 

- in dem die Integrationen vertauschbar sind, ergibt, je nachdem man die innere 
oder dussere Integration zuerst ausfiihrt: 

; ae : 

‘ T(x) = 2 | ud eS ate ges / ge ue Ay 
es ee 

- Sy fo) 

; Of oa eS 

x y2+(k—1) y2(k—1) 

r 0 

Ce 

4 1) Calanda S.A., Domat/Ems; Hispano-Suiza (Suisse) S.A., Genéve. 


2) H. Hésti, Storungen des Raketenfluges durch Windeinwirkung und verwandte Einfliisse, Flug- 


wehr und Technik 1952, 177 ff. k 
3) Das heisst eine relative Genauigkeit von 19/o9 fiir alle vorgesehenen Werte von k. 


— Tr = es Ye eT 
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oder, wenn k = ctgw, das heisst 
e-tw 


1 vs 
ek aaah = Bp 1 
k—i=ctgw—i= =— (cosw —isinw) = = (1) 
gesetzt wird, 
Ps i 1/2 
sin w 
= ug FP ar Foo) ) piele (2) 
V2 x J 
Es ist somit 
x 
Ja) se ft at 
V2% " 
[oe] [oe) 
a ie, (/ e—t(k—i) 4-1/2 ar— f ete? t—l2 dt 
Z. 
ie : L (3) 
[o,0) 
ea a= te te p12 ay 
V2 
x 
[o.e} 
i eae 
= (Se) eter — wats Nala) = fe 1 at 
2 eae | 
x 
2. 
Das Integral 
[@,2) 


N,,(*) = fete pont D2 dt ; 
x 


welches fiir = 0 mit N,(*) tibereinstimmt, geht durch partielle Integration 
(% = tO" DR gy = e—*@—-9 gi) iiber in 


foe) 
1 tea 2 els 2 aye I 1 
N,(*) = |- ag =i) an Ter Le ‘ —t(k—i) ,—(2n+3)/2 
x 
1 i) 2)/2 27 + i 
= 1) eho", en eee 
a ear . 2 ‘aa Nnil*) . 
Es gilt also 
1 —_ —4) 2 1 ul 
N, => — x(k v) —1/2 So 
sli eatarperree a oe Mt), 
Least 2 3 1 
N,(x) = — Sg #1) 8/2 eh 
(7) cl eee a o} Rog Nal4) » 
N,(#) = + Le e-*0-) y-@nsye_  2m+1 1 x) 
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Multipliziert man in diesem Gleichungssystem die erste Zeile mit 


die zweite mit repeat 
as 
die dritte mit 
fed)? As ‘ ‘ : ‘ = usw. 
nd die letzte mit : bebe ; 
(1927: gals m is 


und summiert hierauf samtliche Zeilen, so verschwinden alle Glieder, die N,(x), 
” ..., N,(#) enthalten, und man erhalt: 


Nola) = + e-*@-*) ‘ tye 1 re 2m—1 A 1 yo m+1/2\ 
ey DF ate ae 2 (k — i)m+4 i} 


lee) 


be Bilin 8 2n+1_ 1 (aac p-O2n+ 32 gy 


4 a ey ia (k — i)m+1 
s 


oder mit eine von (1) und der Identitat 


in 


& de alt ao ib 3) 2m—1 2 4 2m (ea ,_1 
2 


. a loin SS | —_., 4 
=4 Aa Z lee a 2 m 22m 4 


n 
Gy = Ey te= 9,5" (0 ica i oa sin™t! wy eimt)wi ean 
m=0 
[e.e) 


=(° ci (n +1)! ew art oy, je Cage j-(2n+3)/2 a4 
n 


x 


Zi jks eer {cd a A =, ae sin” t! w elm+De+ ali aml 


— R,(%) = +Sy(4%) — Ryl*) - 


S,,(¥) wird offensichtlich fiir jedes feste m gleich Null, wenn x gegen Unendlich 
eht. Um zu zeigen, dass unter den gleichen Bedingungen auch das Restglied 
R,,(x) verschwindet, geniigt es, zu beweisen, dass fiir » = 1, 2, 


CO 
lim oe t(h-1) ¢-(2n+3)/2 dp-= 0 
A> Wy 


st, womit dann feststeht, dass S,.(¥) eine asymptotische Reihe darstellt. 
Es gilt nun wegen |e**| =1 und x>1 
a) fir k +0: 


[o.e) 


> en uk 
on ye -C2n+32 qe] — y-(2n+3)/2 a = D,<cund 
d 


er he m/e et we ls be a lh Tad No = Nag Su? ee Sadie th 


mo . 
ie OSes dt\< 


a eh AT ee 
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r 


b) fiir k = 0 ; 
co 

ena? p-(2n+3)/2 gy 3 
2 ; 
Da offenbar sowohl lim L, als auch lim L, bei festem » gleich Null wird, so 


4%—>0O x2 0O F 
konvergiert FR,,(¥) gegen Null, wenn 7 + co geht, was zu beweisen war. 


7 iy oye | 2 2n+1)/2 
<| fi‘ m+2)/2 ay 4 2n+d) =L,. 
x 


| @a+1 


3. 


Die Substitution von (5) in (3) liefert demnach die gesuchte asymptotische 
Darstellung: : 


x 


Rites) = —— | ore) z-1/2 dt ee (ea eiwi2 = 1 e-ke 42 


y2x 2 2x 
(6) 
= 2m 1 ane 
x ee ( = ) m! ak sin™*1 w elim+l)w+a)i xml : 
m=0 


Durch Zerlegung von J(#) in Real- und Imaginarteil ergeben sich die asym- 
ptotischen Entwicklungen 


x 
1 ‘ i 1/2 
ira [coste™™ t~U? dt w cos (=) ate" (2 70 gE | 
7% 
~ (2m 1 
x DS. (—1)™ ( ye ) m! ET sin”*? w cos[(m + 1) w+ #]*~™ 
m=0 
w (sinw)\1/2 (ie (7) 
= cos ( 5 ) = oe a {sinw cos (w + x) 
— sin?w cos(2 w + #) ei + sin’ w cos(3 w + *#) 2 
2H 4 x? 
: ak) P 105 
— sintw cos (4 w + *) ——— 5 ee En 
( + ) ya t+ Sin?w cos (5 w + 4) are + | 
und 
x 
1 ae a] i »\1/2 
——— | sinte-*'-1? dt w sin (= hes Bere: ~ 12 
yaa 5 3 e (2 x0. %) 


CoO 
2m 1 . ¢ 
% aS, (—1)™ ( ss ) m|\ ae sin”*? w sin[(m +1) w+ 4%] «~™ 
m=0 
w (3 


1/2 . 3. 
—'..|—— ) —e** (2 ax) ~ V2 {sinw sin (w + #) 


ie 1 3 
— sin?w sin (2 w + »#) ae + sin®w sin (3 w + x) hice 

a 
105 


; : LS 
— sintw sin (4 w + *) —— + sin5w sin (5 yee 
3 Pata 


8% 
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Setzt man k = 0, das heisst w = /2, so erhadlt man die entsprechenden asym- 
ptotischen Reihen fiir die Fresnelschen Integrale: 


: x 
. 1 1 

C(x) = = | cost #4 dt 5 + 20 x)~ 12 
- x2 aes 


CO 
x Y(-1)" oa me eee Oe 


ep 
_ sche ts 16 x4 


und 


. 
| ag: 


22 


> 


a 1 % + 
Ae D> (-1)" (" ‘x m! cos eg 


* = fosinee? dt = z — (20 x)~¥2 
4 
a 


a 
m=0 si am = (10) 
Ip oe (2 2 x)—1? coae # site — cose : — sin# ee 
eee: 2% 44? 8 x8 
4 se ee ais 
4 + cos Ta : 
E 
j 
; 4. 


, 
B- 


_ Dain asymptotischen Reihen das Restglied bzw. dessen Absolutwert |R,,(4)| 


‘den des letzten Gliedes der Reihe S,,(x) nicht iibersteigen kann, so gilt mit Beach- 
tung von (3), (4), (5) und k = 0, x > 0 die Abschatzung: 


E [Rala)| S feo#* (2 2) 2? (—1) (2) ml ae elim Der al sin wa 
a % 

i Ute aes) 2n—1 

3 2 =1/2 -n 

. <le kx (2 nx) (en ee eee ; x 

; : ee aes ee 

, — 1/2 n 

p < (2 2 x) -e, 3 P ates 3 

: . 

_ Daraus ergibt sich nun fiir 

1 eh ey Paneer ee o 
se *+=10,n=4: Pare 5 a ee de ae ae 20 Eger 
4 Ef 

und fiir 

| 1 Le -2,5 : £4 
- x= 20, n=2: |R,(20)|< tiie cree A pe 


== 
S) 
3 
NS) 
rw) 


.= 7 


ZAMP 


Kurze Mitteilungen — Brief Reports - Communications bréves 


350 


v9 OF FE GFT + O8T 
v9 CT ,LE OTT + 08T 


uS bP ,6E £8 +.08T 
‘uS‘OT,ZF.0S +.08T 
ub 8b ,PhoLT +08T 


un 02 LP ob9T +.08T 
v€ CS Oo TEL + .08T 
ue ¥C 26.86 +.08T 


nL 96 ,46.S9 a o08E *9 


+ 081: 4 = 


ST 


TE OOT 
8é 06 
S~ 08 
co OL 
61 09 
os 
ai OF 
O€ 
0d 


doe = x 


DE Bo rLGacs + ,08T 
u€ EF ,6€ SST +.08T 
uS 85,12 .86 +.08T 
uLET WY olv +.08T° 


NNANM 


16 82,9¥.€9T+ oO8T 4 ji 


v0 bb 82.901 +.08T : T 
ue 66 OT 06h +.08T: T 
uv VL ES oTLE 
19 62,/SE oFLE+ 
ye PRL L ols ot 


»+,081:%=d 


dor =x 


SSDUJAYULY SUL ap4aM-¥ WapuamMmoY advAT ur sap SunuYyrasmpQ 


é Adel 
pee os ee ee ee 
6G L220 SE 6rS‘O 6+ 6¢8'0 10 0 o$22 | 10 ,0 .O8T 
ZO8b7'0 68 966'0 T-9TTL8'0 u€S,E OVE | 19 J€ COL 
890220 6TE9S‘0 T-8S Z68'0 u8 CV o9SS | yoV/TE $07 
LS S670 89L96‘0 T-0r £160 uSE,0 GLE | 9S, 0% 
797ZE‘0 T869S‘0 T-v@€ £6'0 wIT,EL (S62 | wlS,8 .9€2 
8S TSE‘ 8889¢‘0 T-S¢ 1S6'0 u6COToLTE | yEC HE EST 
86 T8E‘0 TZ +9S‘0 T-LLL96‘0 uVE OS OVE | 6€ LY ocLZ 
60€T+r'O €7SCC‘O 1-62 T86'0 wET/OE099E | ,OT,ZT £62 
C6 Eby 0 OS TrS‘O T-8+ 1660 u€ LC o€6E | y8E/SPoFTE 
6ZEL+'0 00€ZS‘0 T-¥8 L66'0 u8¥ 9 olCF | 97,6 oLEE 
00006‘0 0000S‘0 00000‘0 uO ,0 .OSF | 0 ,0 .O9€ 
Gok 
— soo @UIS SOT OG ay 
zit\ @uis} m 


(OT) 


v0 ,O oSET | «O 0 06 v0 ,O0 SF OT 
0,6 obVbT | wEELT 96 | yS‘9F,0 8h 6'0 
w9T,L o¥ST | wTS,0FoZ0T | wS°SZ,0Z oTS 8‘0 
wh0,- oS9T | 85,0 OTT | »0°6Z,0 oSS L‘0 
HEV GLE | 80 o8Ll | no VLC 06S 9'0 
wL1,8To06T | yIT,2S.92T | vL'S ,97.€9 C0 
uV¥,SEobOE | OS ,E€T9ET | w8‘FS,1T 089 +0 
18 VS 0612 uS /9E OFT u9 TC BToEL €0 
uVUb o9ES | vOV/CTOLST | uS'F7,TH BL Z'0 
uS 0S o0SC | uEV/VES89T | ,9'TZ/LT ob8 T‘0 
uO ,0 OLE | 0 ,0 .O8T u9 ,O0 06 0°0 

ag 2S @ cf 


$1q (1) appasaquy sap Sunmmysag UAYISIABUNU ANZ aysamsfjUpT 
T 911992 L 


A 


a Vol. V, 1954 Kurze Mitteilungen — Brief Reports — Communications bréves 351 


-* Die seforderte Genauigheit von + 2-10-4 wird somit bei simtlichen Reihen (7) 
bis (10) erreicht, wenn man fiir 10 < x S 20 bei m= 4 und fiir x>20 bei m = 2 
_ abbricht. 


E _ Fir sehr grosses » gelten gemiss (7) bis (10) die Approximationen 


w ( sin w 


‘a - 
- 1 | 1/2 

is —Rbeplf2 ay * —hx ~1/2 . 

yaa! ot t dt ~ cos 5 5 ) e (2 2 x)~4? sinw cos(w + #) , 
, 0 


x 


S 1 ; Se eal . W /(sinw \1/2 
Bia [sinte he p12 = sin 5 ( : ) e-™ (2 2 x)~ 1? sinw sin (w + x) , 
ny 


t (w = arcctg k); | 


1 ae 
Cx) == ait (22)~*sinx und SiG 4) == is — (22 %)— 1? cosx. 


: Die Nullstellen der Funktion C(%) — 1/2 bzw. S(x#) — 1/2 legen somit fiir 
_ 1 < + in einer beliebig kleinen Umgebung von N a bzw. (N + 1/2) « (N = natiir- 
_ liche Zahl). 

: In analoger Weise nehmen die Funktionen (7) und (8) die konstanten Werte 


5 ae w ee 1/2 reins, os w /sinw \1/2 
i 2 2 ) 2 ea 2 ) 


; fiir sehr grosses x an, wenn dieses in einer beliebig kleinen Umgebung von 
> Nx + (2/2) — w bzw. von N a — w liegt. 


1 5. 
: Die nebenstehenden Tabellen geben die Werte von 
w /sinw \1/2 Fe SOON Ui ers : ; ea 
; cos 5 ( 2 ) und sin (=) HCl oo= OOK SOP seer 10) 
mund die Umrechnung von 7 in Grade fir 7 = 1, 2,..., 10, 20, ..., 100 an. 
Summary 


: Calculating the influence of the wind on the trajectories of rockets (see foot- 
note 2) we find functions of the form 


x 
flay fe ae (i = 1), 
0 

-where, in general, 0< k <1 and # = 10. The present article shows that integrals 
of such a kind can be represented by asymptotic series convenient for the numer- 
ical evaluation with regard to the large values of v. The results are given by 
formulae (6) to (8). Putting k = 0 we get the analogous formulae (9) and (10) 
for FRESNEL’S integrals. To facilitate the computation, two auxiliary tables are 
added to this paper. 


(Eingegangen: 26. Oktober 1953.) 
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Erratum 


Boundary Layer on Rotating Spheroids. By BHASKAR SADASHIV FaDNIs, 
Kharagpur, India (ZAMP 5, Fasc. 2, 156 [1954]}). Equations (2.1), (2.2), and (2. 3) 
obviously contain many more terms on their right-hand side and this should be 
indicated by + --:. In fact, equations (2.1), (2.2), and (2.3) are written down from 
the vector equation 
| ia 
Q 


where v and w are the velocity and the vorticity vectors respectively. 


1 . 
v x w = grad ( +5 04) + ycurl wo, 


Buchbesprechungen — Book Reviews — Notices bibliographiques 


Collected Papers of Peter W. J. Debye. (Interscience Publishers, Inc., 
New York, 1954.) 700 pp., $9.50. 

Es ist einem ehemaligen Schiiler nach Durchsicht der Gesammelten Schriften 
Prof. DEBYES, die anlasslich seines 70. Geburtstages (24. Marz 1954) von seinen 
Freunden und Mitarbeitern herausgegeben wurden, selbstverstandlich eine freu- 
dige Pflicht, die Aufgabe des Rezensenten tibernehmen zu diirfen. Er erinnert 
sich dabei wieder der begeisternden Vorlesungen, welche DEBYE in den zwanziger 
Jahren in Ziirich an der ETH. gehalten hat, ihrer bewundernswiirdigen Klarheit 
und Sachlichkeit und ihrer schépferischen Mittelbarkeit. Was aber vor allem aus 
den Gesammelten Werken des Gefeierten spricht, ist die erstaunliche Vielseitig- 
keit des Gelehrten, der es auf den verschiedensten Gebieten der Physik und 
Chemie verstand, neue Erkenntnisse zu vermitteln. 

Die Gesammelten Schriften sind in fiinf Abschnitte gegliedert. Im ersten be- 
finden sich die Arbeiten aus dem Gebiete der Réntgenstrahlen. Als Kernpunkt 
steht das mit SCHERRER entwickelte Verfahren zur Bestimmung der Kristall- 
struktur aus den Interferenzen von R6ntgenlicht an regellos orientierten Teil- 
chen, welches Verfahren unter dem Namen Debye-Scherrer-Methode fiir die 
Rontgenometrie von ausschlaggebender Bedeutung wurde. Im zweiten Abschnitt 
folgen die Ver6ffentlichungen, die im Zusammenhang mit dem Begriff des elek- 
trischen Momentes der Molekiile stehen. Sie spielen fiir den Chemiker in Fragen 
des strukturellen Aufbaus der Molekiile eine massgebende Rolle. 

Der dritte Abschnitt umfasst die Schriften tiber Elektrolyte, welche gréssten- 
teils in Zusammenarbeit mit andern Autoren entstanden sind. Sie fiihrten vorerst 
zu einer Klarstellung der fundamentalen Begriffe und bilden mit ihren Modell- 
vorstellungen die Basis der modernen Elektrolytchemie. Der vierte Teil bringt 
Arbeiten tiber Lichtstreuung in hochpolymeren Lésungen und kolloidalen Sy- 
stemen. Diese Mitteilungen sind ausschliesslich neuesten Datums; sie gelangten 
aber rasch in der Physik der Hochmolekularen zu anerkannter Bedeutung und 
fihrten unter anderem zu einer einwandfreien Bestimmung des Molekularge- 
wichtes. Im letzten, fiinften Abschnitt sind endlich jene Veréffentlichungen 
DEBYES wiedergegeben, die sich nicht in ein bestimmtes Erscheinungsgebiet ein- 


a 
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_ordnen lassen. Unter ihnen finden wir die Arbeit liber Lichtstreuung durch 

_ Ultraschallwellen, welche die heute unter dem Namen Debye-Sears-Effekt be- 
kannte Erscheinung behandelt. 

____ Das umfangreiche Sammelwerk wird nicht nur den ehemaligen Schiilern und 

_ Mitarbeitern eine grosse Freude bereiten, sondern auch jedem interessierten Leser, 

weil es Einblick in das Leben und Wirken einer Personlichkeit bietet, die es ver- 
standen hat, einerseits in meisterhafter und giitiger Weise den Nachwuchs von 
‘Lehrern und Forschern zu férdern und anderseits ein bleibendes Werk eigener 

_ wegleitender Erkenntnisse zu schaffen. R. Sanger 


Progress in Nuclear Physics, Vol. 2. By O.R. Friscu (Pergamon Press 
Ltd., London, 1952). 295 pp.; £3/3/0. 


Die Kernphysik hat in jiingster Zeit so ungeheure Fortschritte gemacht, 
‘sowohl auf experimentellem wie auch auf theoretischem Gebiet, dass es fiir den 
einzelnen Forscher unméglich geworden ist, alle Teilgebiete in gleicher Weise 
_ zu beherrschen. Deshalb ist auch dieses Buch sehr zu begriissen, in welchem aus- 
gesprochene Fachleute ihre Kenntnisse in einzelnen Artikeln zur Darstellung 
_bringen. Dabei handelt es sich um eine Auswahl besonders aktueller Probleme: 
1. The Electron Optical Properties of Magnetic Beta Ray Spectrometers (N. F. 

VERSTER) ; 

. Nuclear Paramagnetic Resonance (R. V. Pounp) ; 

Luminescent Materials For Scintillation Counters (G. F. J. Gartik) ; 

The Neutron Proton Interaction (G. L. SqurrEs); 

. Fission (W. J. WHITEHOUSE); 

. The Low-Lying Excited States of Light Nuclei (W. E. BURCHAM) ; 

The Nuclear Shell Modell (B. H. Frowers) ; 

. Ionisation by Fast Particles (T. E. CRANSHAW). 

5 Dieses Buch gibt einen ausgezeichneten Uberblick iiber den heutigen Stand der 
Physik auf diesen Gebieten und vermittelt am Ende jeden Kapitels Literatur- 
angaben tiber die einschlagigen Arbeiten des betreffenden Wissensgebietes. Es 
ist jedem Kernphysiker, der sich fiir die theoretische Deutung experimentellen 
Materials interessiert, sicher eine wertvolle Lektiire. O. Huber 
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Nomogramme zum praktischen Rechnen. Von E. V6.LLM, Ziirich. 
Herausgegeben von Ep. ArrniI-LEvcH, Bern. Preis pro Blatt 35 Rp. 

Von den zwei gedruckten Rechenblattern handelt es sich beim ersten um eine 
sehr zweckmassige Tafel zur Bestimmung von logarithmischen Leitern beliebigen 
‘Moduls in der Nomographie. Die zweite Tafel stellt ein Nomogramm fiir den 
Sinussatz der ebenen und der spharischen Trigonometrie dar. 

Das erste Blatt wird jedem willkommen sein, der praktische Nomogramme 
zeichnet, wahrend beim zweiten besonders auf den didaktischen Wert fiir den 
Unterricht hingewiesen sei. A. Kiinzi 


Kernel Functions and Elliptic Differential Equations in Mathematical 
Physics. By S. BERGMAN and M. ScuiFFeR (Academic Press Inc., Publishers, 
New York, 1953). 

In den letzten Jahren hat die von BERMAN und BocHNER 1922 eingefiihrte 
Kernfunktion unter der Zusammenarbeit von BERGMAN und SCHIFFER bei 
Randwertproblemen elliptischer Differentialgleichungen eine zentrale Stellung 
erobert. Ist die zu einer elliptischen Differentialgleichung in einem Gebiet ge- 
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hérige Kernfunktion bekannt, so lassen sich die Lésungen der Randwertaufgab 
erster, zweiter und dritter Art sowie die Greensche, Neumannsche und Robinschi 
Funktion durch Quadraturen herstellen. Der erste Teil des Buches gibt eine Ub 
sicht iiber die bei verschiedenen Gegenstanden (Warmeleitung, Hydrodynam: 
Elektro- und Magnetostatik und Elastizitatstheorie) vorkommenden Randwert 
probleme als Motivierung und Illustration der im zweiten Teil behandelter 
Theorie. Ay. Pfluger 


Aerodynamics of Propulsion. By D. KGcHEMANN and J. WEBER. (McCaul 
Hill Publishing Co. Ltd., London 1953). 331 pp., 150 figs.; $9.-. . S| 

Mit ihrer Publikation fiigen die Verfasser der Schriftenreihe J. C. HUNSAKERR 
«McGraw-Hill Publications in Aeronautical Science» einen weitern Band bei und 
geben damit einem neuen Zweiggebiet der Aerodynamik den Namen. Die bishe- 
rigen Gesamtdarstellungen der Aerodynamik befassen sich vorwiegend mit der 
Erzeugung des Auftriebs, wahrend hier die Str6mungsprobleme der Erzeugung 
des Vortriebs erstmals zusammengefasst werden. 4 

Dem Kernproblem der Vortriebserzeugung, der mechanischen und thermi- 
schen Energiezufuhr an Gasstréme sind die allgemeinen Kapitel gewidmet. We- 
sentlich ist nun aber, dass in die Vortriebsaerodynamik die zahlreichen Probleme 
der Luftfithrungen, Kanalstr6mungen, Einlaufformen, Miindungen, wie sie bei 
Luftentnahmen, Tunnelschrauben, Kiihlern und bei jeder Triebwerksinstallation. 
in Flugzeugen zahlreich auftreten, eingeschlossen werden. Was nur miihsam und 
weit verstreut in der Literatur der Luftfahrttechnik aufgefunden werden kann, ist 
hier in verwertbarer Form und reich belegt durch zum Teil unveréffentlichte 
Messungen geordnet wiedergegeben. Das Buch begrenzt das neue Teilgebiet 
einerseits mit der Strahlausbreitung der Strahlantriebe und andererseits mit den 
Vortriebserzeugung des Insektenfluges. 4 

Die aerodynamischen Probleme der Turbomaschinen und Propeller sind mit 
Recht ausserhalb der Betrachtungen gestellt, und die Propulsionsaerodynamik 
bei Uberschall als noch zu unerforscht nicht in die Arbeit aufgenommen. Eine 
passive Form der Vortriebserzeugung, die Grenzschichtabsaugung, tritt hier im 
‘den Hintergrund, wahrscheinlich weil sie ebensosehr mit der Profilaerodynamik 
wie mit der Vortriebserzeugung verkniipft ist. W. Spillmania 


Proceedings of the Western Computer Conference. Held by the Joint 
ARE-ATEE-ACM Computer Conference Committee Los Angeles, California, Fe- 
bruary 4 to 6, 1953. (Published by The Institute of Radio Engineers, Inc., New 
YOrk,-1953).2231 pp. 135 figs,» $3.50: | 

Die iiberaus schnelle, an vielen verschiedenen Orten betriebene Entwicklung 
der elektronischen Rechengerate macht einen regelmassigen Gedankenaustausch 
der beteiligten Forscher zu einer Notwendigkeit. Der vorliegende Band ist das: 
Protokoll einer solchen Tagung. Digitale und Analogierechengerate sind in gleicher 
Weise beriicksichtigt. Das Hauptgewicht liegt hier nicht so sehr auf den tech- 
nischen Details der Rechengerate als auf organisatorischen Fragen des Einsatzes; 
ferner auf den Erfahrungen bei der Lésung praktischer Probleme, die zumeis+ 
dem Flugzeugbau entstammen. Interessant sind verschiedene Verfahren zu 
Umwandlung zwischen digitalen und Analogiegréssen sowie die standig zuneh} 
mende Verwendung von Netzwerken fiir partielle Differentialgleichungen. Dat 


Buch vermittelt dem Leser einen iiberaus wertvollen Einblick in eine gross¢ 
Menge praktischer Erfahrungen. A. Speise} 
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Mitteilungen aus dem Institut fiir angewandte Mathematik an der- 
Eidgenéssischen Technischen Hochschule Ziirich, herausgegeben von 
Prof. Dr. E. Stiefel. 
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Tabellen zur Erzeugung von Funktionen einer und zweier 
Variablen mit linearen Potentiometern. Von Moheb Aziz 
Abdel-Messih. (1954). 33 Seiten. Broschiert Fr. 4.70. 


An Oscillation Theorem for Algebraic Eigenvalue Problems 
and its Applications. Von Frank W. Sinden. (1954). 57 Seiten. 
Broschiert Fr. 6.25. 


Friiher erschienen: 
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Entwurf eines elektronischen Rechengerates unter besonderer Be- 
riicksichtigung der Erfordernis eines minimalen Materialaufwandes 
bei gegebener mathematischer Leistungsfahigkeit. Von Ambros P. 
Speiser. (1950). 67 Seiten. Broschiert Fr. 6.75. 


Programmgesteuerte digitale Rechengerate (elektronische Rechen- 
maschinen). Von Heinz Rutishauser, Ambros Speiser und Eduard 
Stiefel. (1951). 102 Seiten. Broschiert Fr. 8.85. 


Automatische Rechenplanfertigung bei programmgesteuerten Re- 
chenmaschinen. Von Heinz Rutishauser. (1952). 45 Seiten. 
Broschiert Fr. 5.70. 


A. OSTROWSKI: Vorlesungen iiber Differential- und Integralrechnung. 


Band III: Integralrechnung auf dem Gebiete mehrerer 
Variablen. 475 Seiten mit 36 Figuren. (1954). 
Gebunden Fr. 78.-, broschiert Pr 73.85. 


ne 


# 


VERLAG BIRKHAUSER - BASEL/STUTTGART 


‘ 


SCHUTZEN MATTSTR. = CTE HER : Co BASEL 


